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(i) 


ABSTRACT 


In  this  thesis  we  have  presented  a  theory  for  sums  having 

the  form 

(!)  hM  =  ^f(ax,a2)  f(a2,a5)  ...  f(a._1,a.)  g(a.)  , 

where  the  sum  extends  over  all  compositions  (a^^a^, . . . ,a^)  of  n 

into  an  unrestricted  number  of  positive  parts.  The  numbers  [f(m, n)} 

and  {g(n)}  are  given  by  generating  functions  F(x,y)  and  G(x) 

respectively,  andwe  seek  a  relationship  between  these  functions  and 

the  generating  function  B(x,y)  =  b(a,n)  yaxn  ,  where  b(a,n) 

n,a 

denotes  the  partial  sum  obtained  from  (l)  by  summing  over  just  those 
compositions  of  n  with  a^  =  a  .  We  show  in  Chapter  I  that 

(2)  B(x,y)  =  G(xy)  +  f  F(xy,^)  B(x,s)  ~  ^ 

J  c 

where  c  is  a  contour  in  the  s  plane  which  includes  the  singularities 

of  F(xy,— )/s  ,  but  excludes  those  of  B(x,s)  .  The  relation  (2)  is 
s 

a  special  case  of  the  Fredholm  integral  equation  which  can  be  solved 
by  standard  methods  described  in  Chapter  I  . 

Chapter  II  begins  with  a  concise  definition  of  an  n-omino, 
and  a  description  of  the  work  done  on  the  cell  growth  problem.  Eden  [1] 
proved  that  the  following  special  case  of  (l)  gives  a  lower  bound  for 
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(ii) 


the  number  of  n-ominoes; 

(3)  b(n)  =  (a1+a2-l)(a2+a^-l)  ...  (a^+a^l)  . 

Using  the  methods  described  in  Chapter  I,  we  find  the  generating 
function  of  {b(m))  as  defined  by  (3);  also,  we  find  this  generating 
function  by  an  elementary  method  in  order  to  show,  in  one  instance 
at  least,  the  advantage  of  treating  the  Fredholm  equation. 

In  Chapter  III  we  prove  there  are  more  than  (3.6)n  n-ominoes 
for  all  sufficiently  large  n  .  To  do  this  we  prove  that 

(4)  f(a1,a2,  . .  ,,ai)  >  f(ai,a2)  f(a2,a^)  ...  f(ai_][,ai)  , 

where  (a^,a2, . . . ,a  )  denotes  any  composition  of  n  ,  and  f(p,q,...) 
denotes  the  number  of  n-ominoes  with  p  cells  in  the  first  row,  q  cells 
in  the  second  row,  and  so  on.  Next,  we  find  the  generating  function  of 
{f(m,n)}  ,  and  use  the  Fredholm  equation  to  find  a  lower  bound  for 

(5)  b(n)  =  ^f(al,a2)  f(a2,a^)  ...  f(aj,_1,ai)  , 

which  is  in  turn  a  lower  bound  for  the  number  of  n-ominoes. 

Read  [ 17 ]  stated  and  proved  a  general  combinatorial  theorem 
which  he  used  to  find  the  number  of  n-ominoes  with  1,2, 3,4,  or  5 
rows  of  cells.  In  Chapter  IV  we  solve  Read's  problem  by  means  of  the 
Fredholm  equation.  Also,  we  indicate  a  method  for  finding  the  number 
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(iii) 


of  n-ominoes  with  k  rows  using  linear  difference  equations. 

Harary,  Prins,  and  Tutte  [6]  used  Polya’s  theorem  to  enumerate 
planted  plane  trees;  also,  they  gave  a  one-one  correspondence  between 
the  planted  plane  trees  with  n+2  nodes  and  the  (tri-valent)  planted 
plane  trees  with  2n+2  nodes  in  which  each  node  has  valence  either 
1  or  3  •  In  Chapter  V  we  show  that  the  methods  of  Chapter  I  apply 
to  these  enumeration  problems  as  well,  and  go  on  to  enumerate  special 
classes  of  planted  plane-trees,  for  example,  planted  planes  trees,  having 
nudes  with  valance  either  1  or  k  ,  Furthermore,  we  give  a  very  simple 
correspondence  between  the  planted  plane  trees  and  the  tri-valent  planted 
plane  trees  by  matching  these  trees  with  the  elements  of  a  certain  set 
of  binary  codes . 

We  discuss  the  connection  between  the  Fredholm  equation  and 
the  generating  functions  for  certain  arithmetic  functions  in  Chapter  VI. 
Sums  extended  over  the  divisors  of  n  may  be  thought  of  as  sums  extended 
over  just  those  compositions  of  n  in  which  all  of  the  parts  are  the 
same.  Also,  by  defining 


(6) 


f  (m,  n) 


f*(m,n)  ,  if  m  <  n 

> 

0  ,  otherwise 


the  index  of  summation  in  (l)  is  actually  restricted  to  the  partitions 
(a^  <  a^  <  ...  <  a^)  of  n  ,  since  all  other  compositions  of  n  only 
contribute  0  to  the  sum.  Using  this  idea,  we  are  able  to  establish 
a  link  between  the  Fredholm  equation  and  the  theory  of  partitions. 
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Chapter  ¥11  contains  still  another  application  of  Polya's 

fundamental  theorem,  this  time  to  find  the  number  of  non- isomorphic 

generalized  graphs  with  n  nodes.  A  generalized  graph  is  a  pair 

(Nn,N)  ,  where  =  {l,2,...,n}  is  the  set  of  nodes ,  and  N  is 

some  subset  of  the  k-tuples  N  x  N  X  ...  X  N  called  the  edge  set 

n  n  n  — 2 - 

of  the  graph. 
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CHAPTER  I 


GENERATING  FUNCTIONS 


Our  objective  in  this  chapter  will  be  to  develop  a  theory 
for  treating  certain  sums  which  arise  in  a  variety  of  combinatorial 
problems.  The  notation  and  accompanying  definitions  introduced  at  this 
time  will  be  used  in  subsequent  chapters  without  further  explanation; 
of  course  a  symbol  defined  here  will  have  various  interpretations  later 
on,  depending  on  the  context  in  which  it  is  used. 

Suppose  (f(m,n)  i  m,  n  =  1,  2,  ...  }  and  {g(n)  :  n  =  1,  2,  .. 
are  given  sets  of  numbers  and  consider  the  numbers  b(n)  defined  by 

(1)  b (n)  =  ^  f(ax,a2)  f(a2,a^)  ...  g^)  , 

where  the  sum  extends  over  all  compositions  (a^,a2, . . . ,a^)  of  n 
into  an  unrestricted  number  of  positive  parts.  We  make  the  convention 

i 

that  g(n)  is  the  contribution  to  the  sum  when  the  number  of  parts  of 
the  composition  is  one. 

The  symbol  b^  (a,n),  used  with  all  or  only  some  of  the  suffixe 
will  denote  the  partial  sum  obtained  from  (l)  when  the  index  of  summation 
has  been  restricted  to  those  compositions  of  n  which  have  exactly  k 


o. 


»  qalsvab  03  •»  J.'.H  '  1  >[,d©  *u0 

. 

:  ,  JJ  13  :  :3  IU':.  i.  <  .-.US'.  >!•>*<  '*  ’’  '  1  ’  ‘  ' 

■ 

M  •  i  a  :  (*!'  :'  fc"»®  (  •  t  t-  *  f  5  0 v  f  '' iC'l'q'  " 

Xd  b^»b  («)<r  8i»4Ptun  3rii  3»biano3  bn»  8i»<teon  io  t3»8  fiavig  •« 


a  lo  (  MoWiac-poa  H»  l«vo  abn  J*o  >«  •*  •>»<* 

1 

io  «3iBq  lo  lodmin  a<U  •>•>!*  *«»  »r  I  ol  noijudlwnoo  »«:’->  •!  M*  J“ 


a0M  ilu*  ’  -on.  ,  J  Ofc-Sj  (I)  40*1  bo-  Ido  «u«  •*»  9,0,  sb  lr" 

>  XI3  <:xs  ov  d  ilsi/lw  «  1  i  'ioi :  j»r,  »»o(U  »3 


-  2  - 


parts,  no  part  greater  than  j  ,  and  the  first  part  equal  to  a  \  if 
a  suffix  is  dropped,  the  corresponding  restriction  on  the  index  of  summation 
is  dropped  as  well.  This  definition  implies  certain  obvious  relationships 


between  the 

numbers  b^(a,n)  defined. 

For  example,  we  have 

n 

n 

(2) 

b(n)  *  \  b(a,n) 

=  Y  bk(n)  , 

L _ 1 

a=l 

k=l 

(3) 

(4) 


(a,n) 


i 


Less  trivial  perhaps  are  the  following  relationships  which 
play  a  central  role  in  later  work:  for  a  <  n  , 


f(a^^^,a_^)  g(a^) 


s*  ^  f(a,v)  b(v,n-a)  , 
v 

(6)  bfc+1(a,n)  =  V  f(a,v)  bk(v,n-a)  ,  k  =  1,  2,  ...  . 

V 

The  index  of  summation  in  the  first  sum  in  (5)  is  extended  over  all 
compositions  (a^a  ,  . . .  ,at)  of  n-a  into  an  unrestricted  number  of 
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positive  parts;  the  second  equality  follows  on  summing  over  compositions 
of  n-a  with  a1  =  v  for  v  =  1,  2,  . ..,  n-a  .  The  proof  of  (6)  is 
similar . 

Now  we  turn  our  attention  to  finding  the  relationships  between 
the  generating  functions  of  the  sequences  being  studied.  Thus,  we  define 

00  00  00 

G(x)  a  g(n)  Xn  , 
m=l  n=l  n=l 

oo  n 

^  b (a,n)  yaxn  , 
n=l  a=l 

so  that  F(x,y)  is  an  analytic  function  of  x  for  fixed  y  and  of 
y  for  fixed  x  in  neighborhoods  of  x  =  0  and  y  =  0  respectively, 
and  G(x)  is  an  analytic  function  in  a  neighborhood  of  x  =  0  .  For 
fixed  natural  numbers  j  and  k  we  write  B^(x,y)  and  B^(x,y)  for 
the  generating  functions  of  {b"^(a,n)J  and  (b^(a,n)}  respectively. 

For  example,  using  this  notation,  the  relations  in  (2),  (3)  and  (4) 
imply 

00 

B (x, l)  =  ^  b(n)  xn  , 

ifci 


(7) 


F(x,y)  =  V  V  f(m,n) 


m  n 
x  y 


(9) 


Bk(x>y) 
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All  of  the  B  functions  are  analytic  if  F  and  G  are 
analytic;  to  show  this  it  is  sufficient  to  show  that 

CO 

(10)  X  X  lf^al,a2^l  lf^a2,a3^l  •••  lf(ai-l>ai)  lls^j.)  I  ^ 

n=l 

converges  for  sufficiently  small  |x|  ,  where  the  inner  sum  extends  over 
all  compositions  of  n  .  The  proof  follows  from  the  fact  that  series 
such  as  in  (7)  converge  if,  and  only  if,  there  exist  positive  constants 
a,  3  and  7  such  that  |f(m,n)|  <  <xm|3n  and  |g(n)|  <  yn  ;  assuming 
this  we  have  for  a  given  composition  (a^a^, .  • .  ,a^)  of  n  , 

(11)  |f(«1,a2)||f(a2,a3)|  ...  |f (ai^1,ai)  |  |f (ai) |  <  (a37)n  . 

n.  ^ 

There  are  2  compositions  of  n,  so  the  coefficient  of  x  in 

(10)  is  less  than  (2a|37)n  ;  hence,  the  series  converges  for  |x|  <  l/(2aP7)  . 

It  is  clear  that  the  B  functions  are  uniquely  defined  by  the 
coefficients  b^(a,n)  once  F  and  G  are  specified;  however,  there 
are  other  ways  to  construct  these  functions.  We  begin  by  finding  a 
relationship  between  F,  G  and  the  sequence  {B^}  • 

By  definition 

r  g(n),  if  a  =  n 

(12)  bx(a,n)  =  4  , 

^  0  ,  otherwise 


in  terms  of  generating  functions  this  becomes. 
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(!3)  B1(x,y)  =  ^  g(n)  xnyn  =  G(xy) 

n=l 

Next,  multiplying  through  (6)  with  yax°  and  summing  over  a 
and  n  we  obtain 

f11*)  Bk+l(X,y)  =  Y_,  X  f^a,V^  bk^v,n"a^  yfl ^ 

a,n  v 

It  is  easy  to  verify  by  substituting  the  power  series  involved  that  the 
sum  in  the  right  member  of  (14)  is  the  residue  with  respect  to  s  of 
F(xy,“)  Bk(x,s)/s  |  hence,  for  k  =  1,  2,  ...  , 

(15)  Bk+l(x'y)  =  f  F(xy's)  VX’s)  T  ’ 

C 

where  c  is  a  contour  in  the  s  plane  which  includes  the  singularities 
of  F(xy,~ )/s  ,  but  excludes  those  of  B,  (x, s)  •  In  practice  the 
singularities  in  the  s  plane  of  F(xy,j)/s  may  be  known,  in  which  case 
the  integration  in  (15)  may  be  carried  out  to  obtain  a  recurrence  relation 
satisfied  by  the  sequence  {B^}  .  If  we  are  able  to  surmise  the  form 
of  B^  in  general  after  calculating  B^,  B^,  ...  ,  this  recurrence 
relation  can  be  used  to  prove  the  conjectured  formula  by  induction. 

Summing  on  k  in  (15)  we  obtain 

00  00 

<l6>  X  Vl(x’y)  =  /  F(xy’s)  I  Bk(x,s)  T  ; 

k=l  c  k=l 

B^(x,y)  =  G(xy)  to  each  side  of  the  equality  in  (16)  we  obtain 


adding 
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(IT)  B(x,y)  *  G(xy)  +  jf  F(xy,j)  B(x,s)  ~~ 

c 

It  is  also  possible  to  find  (17)  directly  from  (5);  here  the  proof 
involves  verifying  that 


(18) 


5a  f  Ifor.jr)  »(x-s)  T 

c 


«>  n-1 


I  I 

n=2  a=l 


n~a 

V  f(a,v)  b(v,n-a)  yV1 
v-1 


y 


but  this  simply  amounts  to  finding  the  coefficient  of  s  ^  in  the 
Laurent  expansion  of  F(xy,— )  B(x,s)/s  „ 

S 


Just  as  in  (15)#  the  integral  equation  in  (17)  becomes  a 
functional  equation  satisfied  by  B(x,y)  when  the  singularities  of 
F(xy,—)/s  are  known;  sometimes  it  will  be  possible  to  use  this 
relationship  to  solve  for  B(x,y)  explicitly. 


i 

Now  we  are  ready  to  consider  the  functions  BJ(x>y)  .  First 

0 

we  note  that  b^(a^n)  would  correspond  to  the  number  b(a,n)  which  would 
result  if 


(19) 


Fj(x,y) 


j 


m=l 


j 

y  f(m,n)  xmyn 
n»l 


and 


Gj(x) 


were  used  in  place  of  F  and  G  to  define  the  sets  (f(m,n))  and 
fg(n))  •  This  observation  along  with  (17)  implies  that 


(20) 


B^(x,y) 


G^(xy)  + 


1 

2iti 


c 


y 


. 
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where  c  includes  the  origin  (since  F(xy,“)/s  is  simply  a  polynomial 
in  1/s  )  but  excludes  the  singularities  of  B^(x, s)  .  Substituting 

•j 

the  representation  of  F„(xy,-)  given  by  (19)  into  (20  )  we  obtain 

J  s 

after  integrating 


(21) 


(x#y)  =  ^  -jg(m)  +  f(m,n) 

m=l  n=l 


nl 


m  m 
x  y 


where  aV  denotes  the  n^  partial  derivative  with  respect  to  s 
of  B^(x,s)  at  s  3  0  «  But  we  also  have 


(22) 


Bj(x,y) 


&V  m 
■  y 


m; 


m=l 


so  that  equating  coefficients  of  yP  in  (21 )  and  (22)  gives  the  system 
of  equations 


(23) 


)PBj 

P' 


j 


g(p)  xP  +  ^  f(p,n)  ^7 
n-1 


x 


for  p  a  1,  2,  j  .  The  system  in  (23)  consists  of  j  equations, 

linear  with  respect  to  the  functions  dPB^/p!  ,  so  Cramer’s  rule  may 
be  applied  to  solve  for  these  functions  explicitly  in  terms  of  x  and 
the  numbers  g(n)  and  f(m,n),  m,n  =  1,  2,  j  .  In  fact,  we  have 


(2*0 


— f--1  =  A  (x)/A  (x) 

Pi  jPv  '  V 


(25) 


A^(x)  =  det  [arg]  ,  Ajp(x)  =  det  [bPg]  , 


■ 


' 
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where  [a  ]  is  defined  by 
rs 


r  f(r,r)  x£  -  1  ,  if  r  =  s 

(26)  ars  =  J 

4  f(r,s)  xr  ,  if  r  4  s 

and  [bPg]  is  obtained  from  [arg3  by  replacing  the  column  of 

[ar£J]  with  the  column  vector  [-g(l)x,  -g(2)x2,  ...,  -g(j)x^]  . 
Substituting  the  expressions  for  8PB^/p*  given  by  (24)  into  (22) 
gives 

j 

(27)  -  zrl^j  £  %(x)  yT"  ’ 

m=l 

4 

so  that  evidently  BJ(x,y)  is  a  rational  function. 


Part  of  our  interest  in  the  functions  b^(a,n)  is  evinced 
by  the  fact  that  bJ(a,n)  =  b(a,n)  for  n  <  j  ,  and  hence  the  implication 
that  B^(x,y)  tends  uniformly  to  B(x,y)  in  a  suitable  domain  D(x,y); 
this  follows  from  the  fact  that 


oo  n 

(28)  B(x,y)  -  ®^(x,y)  »  ^  (b(a,n)  -  b^(a,n))  yaxn  . 

n=j+l  a=l 


In  some  of  the  problems  we  study  later,  we  will  have 
0  <  b^(a,n)  <  b(a,n)  with  (fo^(a,n)  :  n  =  1,  2,  ...  )  an  increasing 
sequence  for  each  j  .  For  example,  this  will  be  the  case  when 
f(m,n),  g(n)  >  1  ,  for  m,n  =  1,  2, 


o  o  o 
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If  P  and  Q  are  polynomials  such  that  Q/P  is  irreducible 
and  generates  an  increasing  sequence  (a  }  ,  then  lim  (a  )^/n 


n  — >  oo 


exists  and  is  the  largest  real  root  of  P(l/x)  =  0  ;  a  proof  of  this 
may  be  found  in  Salem  [21]  .  Thus,  if  B^(x,l)  =  L  AJp(x)/A^(x) 
generates  an  increasing  sequence  (b^(n)  l  n  =  1,  2,  „ . ,  ]  with 
0  <  b^(n)  <  b(n)  ,  then 


(29) 


(b(n))^81  >  (b^(n))^n  >  0 


for  all  sufficiently  large  n  ,  where  0^  is  the  largest  real  root 


of  Aj (l/x)  =  0  .  If  furthermore,  (b^(n)  :  j  =  1,  2,  ...  }  is 


an 


increasing  sequence  for  each  n  (this  will  also  be  the  case  if 
f(m,n),  g(n)  >  1  ,  m,n  =  1,  2,  ...  )  then  0^  <  6^  <  .  <,  <>  is  bounded 
by  lim  (b(n))^/n  =  p  so  that  (0  }  tends  to  a  limit  0  <  p  . 

n  — *  oo  J 


n  oq 

Thus,  for  every  e  >  0  we  have 


(30) 


(b(n))Vn 


>  0 


for  all  sufficiently  large  n  .  So  far  we  have  been  unable  to  discover 
whether  0  =  ft  under  these  assumptions. 

The  relation  in  (17)  is  a  special  case  of  an  important  integral 
equation  known  as  the  Fredholm  equation;  Whittaker  and  Watson  [24]  give 
the  following  definition.  Let  to(z)  be  a  continuous  function  and  suppose 
K(z,s)  is  a  real  function  of  z  and  s  such  that  either  (i)  it  is 
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a  continuous  function  of  both  variables  in  the  range  a  <  s,  z  <  b  , 
or  (ii)  it  is  a  continuous  function  of  both  variables  in  the  range 
a  <  s  <  z  <  b  and  K(z, s )  *  0  when  s  >  z  .  Now 

n  b 

(31)  0(z)  =  ^(z)  +  /  K(z,s )  J2f(s)  ds 

a 

is  Fredholm's  equation;  K(z,s)  is  called  the  kernel.  To  see  that  the 
equation  in  (I7)  has  this  form,  we  suppose  c  is  defined  by  s(t) 
for  0  <  t  <  1  ,  regard  x  as  a  fixed  number,  write  y  =  s(z)  , 
and  convert  the  contour  integral  into  a  line  integral  to  obtain 

(32)  B(x, s (z) )  =  G(xs(z))  +  “t  J*  F(xs(z),^j)  B(x,s(t))  . 

o 

This  equation  is  less  general  than  the  Fredholm  equation  in  (31)  principally 
in  that  F  and  €  are  assumed  to  be  analytic;  consequently,  the  theory 
of  the  equation  is  considerably  simpler  in  the  case  we  consider.  For 
example,  the  existence  of  a  unique  solution  is  guaranteed  since  the 
definition  of  b(a,n)  in  (5)  and  the  relationship  between  the  generating 
series  in  (17)  are  equivalent  statements. 

The  expression  for  B(x,y)  given  by  (9)  corresponds  to  the 
solution  of  (32)  given  by  the  method  of  successive  approximations  and 
is  called  the  Neumann  series  while  the  approximations  of  B(x,y)  given 
by  the  functions  BJ(x,y)  correspond  to  the  solution  given  by  Fredholm's 
method.  An  excellent  exposition  of  the  theory  of  this  integral  equation 
as  well  as  a  description  of  the  work  done  since  Fredholm  first  gave  a 
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complete  theory  for  them  in  1900,  may  be  found  in  Riesz  and  Sz.-Nagy  [ig]. 

The  connection  between  the  Fredholm  equation  and  convolutions 
of  the  kind  given  in  (5)  seems  not  to  have  been  observed.  As  will  be 
seen  in  subsequent  chapters,  certain  arithmetical  functions  as  well  as 
combinational  functions  satisfy  convolutions  of  this  type. 


r  .  '  4fc.M  i  i  a-  VM  lOCK  ■  ■"'* 

,  >•  •  l***W*.  •*»  fl93V**d  flOl^SflflOO  •** 

„  ,  .  t  ,:  r  3-:  Jl  «3  -  <£.  .313*  '.  *.**»  ,«•  <!»  ’  “•"P**4  ; 


12  “ 


CHAPTER  II 


CELL  GROWTH  PROBLEMS 


The  square  lattice  in  the  plane  is  defined  by  Hilbert  [10] 
as  the  set  of  all  points  of  the  plane  whose  Cartesian  coordinates  are 
integers,  A  cell  of  the  square  lattice  is  a  point  set  consisting  of 
the  boundary  and  interior  points  of  a  unit  square  having  its  vertices 
at  lattice  points.  An  n°omino  is  a  union  of  n  cells  which  is 
connected  and  has  no  finite  cut  set . 

The  set  of  all  n~ominoess  R  ,  is  an  infinite  set  for  each 
n  J  however,  we  are  interested  in  the  elements  of  two  finite  sets  of 
equivalence  classes,  and  ,  which  are  defined  on  the  elements 

of  R^  as  follows?  Two  elements  of  R^  belong  to  the  same  equivalence 
class  (i)  in  S  ,  or  (ii)  in  T  ,  if  one  can  be  transformed 
into  the  other  by  (i)  a  translation,  or  (ii)  by  a  translation, 

rotation,  and  reflection  of  the  plane.  An  element  of  R  is  in 
standard  position  if  it  is  above  the  x-axis  with  a  cell  at  the  origin 
and  all  cells  in  the  first  row  are  to  the  right  of  the  y-axis. 

There  is  exactly  one  element  in  each  equivalence  class  in 

S  which  is  in  standard  position,  while  an  equivalence  class  in  T 
n  n 

may  contain  as  many  as  eight  elements  all  in  standard  position.  Thus, 
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if  s(n)  and  t(n)  denote  the  number  of  elements  in  S  and  T 


n 


n 


then 


(1) 


g  s(n)  <  t(n)  <  s(n) 


Harary  [ 4  ]  has  listed  the  cell  growth  problem  as  an  unsolved 
problem  in  the  enumeration  of  graphs.  Stated  in  the  terms  we  have  just 
defined,  Harary' s  formulation  of  the  cell  growth  problem  is  to  find 
t*(n)  ,  the  number  of  equivalence  classes  of  T  which  contain  simply 
connected  n~ominoes,  If  an  n-omino  is  not  simply  connected,  in  the 
sense  that  it  has  "holes",  then  it  is  said  to  be  multiply  connected, 
Harary  [ 5  ]  later  reported  that  a  computer  had  been  programmed  to  find 
t(n)  for  n  <  12  ,  and  he  listed  these  values.  Evidently  this  work 
was  carried  out  independently  by  Stein,  Walden,  and  Williamson  at  the 
Los  Alamos  Scientific  Laboratories  and  by  Lander  and  Parkin  at  the 
Aerospace  Corporation.  This  and  more  detailed  information  about  the 
counts  were  communicated  to  the  author  by  Mr.  Parkin. 

Read  [ly  ]  calls  representative  elements  of  the  equivalence 
classes  of  S^  and  fixed  and  free  animals  with  n  cells 

respectively^  also,  he  gave  a  method  for  finding  the  number  of 
equivalence  classes  of  S^  which  contain  n-ominoes  in  standard  position 
having  cells  in  exactly  k  rows  above  the  x-axis.  He  calculated  these 
numbers  for  k  <  5  and  n  <  10  and  used  these  results  to  find  t(n) 
and  t*(n)  for  n  <  10  |  an  error  in  his  calculations  involving  t(l0) 
was  discovered  by  Parkin  and  subsequently  corrected  by  Read,  The  known 


I 
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values  of  t(n)  and  t*(n)  are  as  follows: 


n 

1 

2 

3  4 

5 

6 

7 

8 

9 

10 

11  12 

(2) 

t  (n) 

1 

1 

2  5 

12 

35 

108 

369 

1285 

4655 

17073  636OO 

t*(n) 

1 

1 

2  5 

12 

35 

107 

363 

1248 

4271 

• — — — •  -  ■■ 

Eden 

[1  ] 

seems 

to 

have 

been 

the 

first 

person 

to  give  upper 

and 

lower  bounds 

for 

t  (n) 

0 

9 

his 

bounds 

are 

(3) 

(3.1^)” 

1  <  t( 

\  ^  1  n 

n)  <4 

9 

for  sufficiently  large  n  ;  his  argument  for  the  upper  bound  is  in  doubt. 

We  will  show  that  a  =  lim  (©(n))^31  exists,  so  that  from 

n  — » oo 

(l)  we  can  conclude  that  lim  (t(n))1//n  exists  and  is  equal  to  a  ; 

n  ->  co 

to  do  this,  we  will  require  a  lemma  due  to  Fekete  (see  Pdlya  and  Szego 
[15],  page  171,  for  similar  results). 

Lemma  1 :  If  *-s  a  sequence  of  natural  numbers  such  that 

f(U  )Vn]  £s  bounded  and  U  U  <  U  ,  then  lim  (U  )Vn  exists, 
n  m  n  ~  m+n  n 

n  -*  00 

To  show  that  (s(n)}  satisfies  the  conditions  of  Lemma  1 
we  prove  two  more  lemmas  £  in  the  proof  of  Lemma  3  we  indicate  the  method 
Eden  used  to  establish  the  upper  bound  given  in  (3)0 

Lemma  2:  s  (m.)  s(n)  <  s(m+n)  ,  m,n  =  1,  2,  ... 


lawoHoi  aft  baa  (n)$  Jo  etulav 


i  oocioq  q.:T  3t  ±*fte  nstd  svari  o}  8ras*» 

. 

'•  'i  V  '  •  ■ 

ai  b  a  eiaixe  n^{(a)i)  mil 


c,^r  >2  .  i*  ay  • ■  i  jv)  $  -'£>:&ri  <  J  !-  fir. a&l  a  it  JIlv  ew  .ttX/l*  ob  03 

.  I  I  '£  -■  •  >“  , 


ir  I  i  jp  ’  Xo  !  o  3,  vat  a  1 


lAV» 
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Proof l 


Let  X  and  Y  be  representative  elements  from  equivalence 


classes  in  S  and  S  respectively,  such  that  the  lower  edge  of  the 

first  cell  in  the  bottom  row  of  Y  is  joined  to  the  upper  edge  of  the 

last  cell  in  the  top  row  of  X  .  The  (m+n)“Qmino  just  described  is  a 

representative  element  of  an  equivalence  class  in  S  *  The  existence 

m+n 

of  this  one-one  correspondence  between  S  X  S  and  a  subset  of  S 

m  n  m+n 

imples  the  lemma,, 


Lemma  3s  s(n)  <  (2j/b)n 


Proof: 


Following  Eden,  we  assign  a  unique  sequence  of  binary  digits 


to  each  element  X  of  which  will  be  denoted  by  W(x)  »  {(cx^p^) 

(apP_70)  . „ ,  (a  P  7  ) }  .  To  do  this,  we  assume  X  is  in  standard 
position  and  draw  a  labelled,  directed  tree  over  the  cells  of  X  ; 
the  nodes  of  the  tree  will  be  the  centers  of  the  cells  of  X  and  the 
directed  edges  will  be  drawn  between  some  of  the  nodes  in  connected 
cells o  The  center  of  the  cell  located  at  the  origin  is  given  the  label 
,  and  directed  edges  are  drawn  to  the  centers  of  the  cells  connected 
to  it l  these  nodes  are  called  C^,  ».*,  proceeding  clockwise  from  the 
y-axis  around  *  Now  we  describe  a  process  that  must  be  carried  out 

at  the  nodes  C^,  C^,  Cr,  the  nodes  being  taken  in  this  order* 

Suppose  the  process  has  been  carried  out  at  C^,  C .*„,  Ch  and 

that  a  directed  edge  has  been  drawn  from  to  „  Beginning  with 


the  side  of  C.  to  the  right  of  C.  and  proceeding  clockwise  around 

i  J 

C±  we  draw  a  directed  edge  to  any  cell  connected  to  (L  which  has 


.  I  aiiO-j  ‘«G  aifJ-  1 


"(ATS)  >  <«).  ■  *5_S2SJ 


„  ,  W  K  .  J>  ii*>  «  lo  «T.JO»»  •*»  *<»  Hiw  ..*1  .&  Jo  »  ‘  or.  «* 

9rf3  os  BWB-rb  «B  I>»sb»  b  ooiib  «  ,  r> 


-  16  - 


not  already  been  labelled,  giving  the  labels  C,  ,  C  1?  to  the 

k  k-fl 

new  nodes,  where  k  is  the  smallest  index  not  used  previously,,  Now 
ai>  ^ i 9  an<*  are  assi§ned  the  binary  digits  1  or  0  if  an  edge 

going  from  cuts  the  first,  second, ot  third  side  encountered, 

proceeding  clockwise  from  around  . 

The  units  which  appear  as  digits  in  W(x)  correspond  to  the 

edges  of  the  tree  drawn  over  the  cells  of  X  ;  thus,  the  binary  digits 

of  W(x)  must  sum  to  n-1  ,  since  a  tree  with  n  nodes  has  exactly 

n-1  edges.  The  number  of  binary  sequences  of  length  Jn-l  that 

3n-l 

contain  exactly  n-l  ones 


/  3n-l  \ 

is  f  o  since  different  elements 

\  n-1/ 


of  give  rise  to  different  binary  sequences,  it  follows  that 


oo 


/  5n*1  \ 
w  5  C ... ) 


*  \  n 

<  (*l\ 


\  *  J 


This  completes  the  proof  of  Lemma  3  J  in  Figure  1  we  show  a 
7-omino  in  standard  position  with  its  labelled  tree  and  corresponding 
binary  sequence. 


(ii )  (on)  (ooo)  (ioo)  (oio)  (ooo)  (ooo) 


Figure  1 


-  • 


»  ol  b  oM.-»noc  (X>  nl  eilgib  .>  9 


b  «...  J-0  •«»  .:-  r  Kit.  *1  1  »*»  '■  1  '  ■  *  0 


■ 


- 


Combining  Lemmas  1,  2,  and  3  we  have  the  following  result. 


Theorem  1.  lim  (s(n))l/n  =  a  exists,  and  a  <  27/4 


n  —>  oo 


If  there  is  a  directed  edge  from  C9  to  C.  .  and  from  C. 

i  j  J 

to  ,  then  the  choice  of  the  binary  digits  in  aac* 

(a.3.7  )  usually  imposes  restraints  on  the  choice  of  (a,  (3  y  )  . 
j  J  j  k  k  k 

Using  this  idea,  Eden  claims  to  have  shown  that  fewer  than  4°  of  the 
binary  sequences  of  length  3n~l  have  digits  which  satisfy  these 
restraints „ 


.th 


If  an  n-omino  in  standard  position  has  exactly  a^  cells  in 


the  j”*‘  row  above  the  x~axis,  for  j  =  1,  2,  ...,  i,  we  say  the 
equivalence  class  it  represents  in  S  belongs  to  the  set  F 


n 


ala2°°°ai 


Eden  observed  that  if  f (a^,a^> » . .,a^)  denotes  the  number  of  elements 


in  F  , 

fl ^ o  O 


then 


(5) 


f(a1>a2»--->a1)  >  (a1+a2-l)(a2+aj-1)  •••  (ai.1+ai*1) 


The  product  in  the  right  member  of  (5)  gives  the  number  of 


equivalence  classes  in  F 


which  contain  n°ominoes  that  have 


12' 


a  connected  strip  of  a^  cells  in  the  j 


th 


row,  for  j  =  1,  2, 
when  the  n-omino  is  in  standard  position.  This  follows  from  the  fact 


that  a  strip  of  r  cells  can  be  joined  above  a  strip  of  s  cells  in 

r+s-1  ways.  Since  there  is  a  set  F  corresponding  to  every 

a  ...ai 

composition  of  n  ,  we  have 


fct  MS-  •  .*  •*  ■  "  *  8Blfli<faa0) 


,  t  3  03  .  W>T*  nbo  bMMllb  •  •*»r:3  11 


...  t:  DT  bW)«  nl  oriiJTo-n  i  M 


bWjt 


,  ,  ,  ,  J  t.  .  ,3  ,J  ,  t  *ol  9rf3  «o*  t  »■=» 

,,  »ri  o  ,  i  i  "-1  s3n»a«q»t  31  se.Io  *M»X»  0p< 


?  t  ore.  ;C  .oolliaoq  amba.31  nl  •»  onlao-n  »«  l  • 
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(6) 


s(n)  >  b (n)  =  ^  (a1+a2-l)(a2+a^-l)  ...  (aA 


•,  +a .  ■ 

■1  l 


•i) 


where  the  sum  extends  over  all  compositions  of  n  into  an  unrestricted 
number  of  positive  parts. 

Eden  was  able  to  show  that  for  sufficiently  large  n  , 
b(n)  >  (3.1^)n  and  in  this  way  proved  the  lower  bound  given  in  (3). 
Klarner  [12]  later  improved  this  estimate  to  (3.21)°  >  b(n)  >  (3.20)n  , 
for  sufficiently  large  n  .  We  will  now  use  the  method  outlined  in 
Chapter  1  to  obtain  the  generating  function  of  [b(n)j  ,  and  at  the 
same  time  prepare  the  machinery  for  treating  similar  problems. 


The  sum  in  (6)  has  the  form  of  (l.l)  with  g(n)  =  1  and 
f(m,  n)  =  m+n-1  J  in  fact,  we  can  suppose  the  numbers  f(m,n)  have 
the  form  w(m)  +  t(n)  )  thus,  if  W  and  T  generate  (w(n)}  and 
{t(n))  respectively,  we  have 


F(x,y) 


m  n 


(w(m)  +  t (n) )  x'y 


yw(*)  +  xrfy) 

1-y  l~x 


In  general,  if  the  kernel  of  the  Fredholm  equation  in  (l.3l) 
has  the  special  form 


(8) 


K(x,y)  =  )  Pk(x)  ?k(y)  , 


b»*a  jr  iin  n  oaol  a  •  >  wolll*o*o>o  II«  two  eboolxs  <m>«  oils  »*«* 


arctics  »  >svo*c  ?'*«•* 


t  11  w  »W  .  a  sh  '•  ^  *  •• 


tu  1  {(o>»)  93  .t*n»3  '  1-0  W  51  »•““«  i  +  '■)«  °*<>-  ' 

■  ‘  .  ••■  « 
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it  is  said  to  be  of  finite  rank  and  the  study  of  the  integral  equation 

reduces  to  a  system  of  linear  algebraic  equations.  (For  more  details 

on  this  point  see  Riesz  and  Sz-Nagy  [18]  ,  page  161.)  Substituting  the 

expression  for  F(xy,— )  given  by  (7)  into  (I.I7)  we  obtain  an  integral 

s 

equation  having  a  kernel  with  finite  rank: 


(9) 


B(x,y)  =  G(xy)  +  W(xy)  B(x,l)  + 


2jti 


I 


xyT (~)  B(x,s)  ds 

b 

(l-xy)s 


Multiplying  equation  (9)  by  T(l/y)/y  and  integrating  we 
obtain  a  relation  which  implies 


do) 


J_  r  Tfi)  b(x  v)  ^  -  p(x)  +  B(*,l),ft,(x) 

2x1  J  T{y>  y  ~  1  -  T(x)  > 


where 


00 


(u)  =  jjJJ  f  ^  =  y  *(”)  s(n)  xn  , 

c  n= 1 


CO 


(12)  Q(x)  =  “  |  T(i)  W(xy)  ^  ^  t(n )  w(n)  x° 

c  rfci 


The  integral  representation  of  sums  such  as  appear  in  (ll) 
and  (12)  was  probably  first  discovered  by  Hadamard  (see  for  example, 
Titchmarsh  [22],  pages  157-159).  Substituting  the  expression  for  the 
integral  given  by  (10)  into  (9)  and  setting  y  =  1,  we  obtain  a  linear 
equation  in  B(x,l)  ;  solving  this  we  obtain 


...  »• 

k  jiKBSiqx* 


15  <••»«  '  *  3 "*’•’» r  noilbUf 


■ 

,  I  *  (-•  r)T  (7  «*M*»P*  >»*  1  * 


.. 


naaiiJ  »■■  -  X-  "  1  '  (  ■’>  i",“ 

,  «  ,,,  ,  a  q  !i  -  jad.*  * s88<!  JSSJ  '  ' 


(13) 


B(x, 1) 


l-x)  G(x)  (l-T(x) 


The  relations  in  (9),  ( 10 ) ,  and  (13)  can  be  combined  to  find 
B(x,y)  in  closed  form  in  terms  of  W,  T,  G,  P,  and  Q  .  When 
G(x)  =  x/(l-x),  then  P(x)  =  T(x)  ,  and  (13)  reduces  to 

(1*0  B(x,l)  =  x/((l-x)(l-T(x))(l-W(x))  -  xQ(x) ) 


An  elementary  proof  of  (14)  can  be  given  as  follows?  Substituting 
f(m,n)  =  w(m)  +  t(n)  into  (1.5),  we  have  for  a  <  n  , 


(15) 


b (a,n) 


w(a) 


n-a 


I 

v=l 


b(v,n-a)  + 


n-a 

^  t(v)  b(v,n-a) 
V=sl 


n-a 

*  w(a)  b(n-a)  + 

v=l 


t(v)  b(v,n-a) 


o 


Writing  a-k  and  n-k  in  place  of  a  and  n  in  (15)  we 
obtain  a  similar  expression  for  b(a-k,n-k)  ;  taking  the  difference 
b(a,n)  -  b(a-k,n-k)  and  transposing  a  term  gives 


(16)  b(a,n)  =  b(a-k,n-k)  +  [w(a)  -  w(a-k)]  b(n-a)  . 

When  k  =  a-1  in  (16),  we  find  that  each  of  the  numbers  b(a,n)  can 
be  written  in  terms  of  b(l,v)  and  b(v)  ,  v  =  1,  2,  ...  \  thus, 


t(x-l)\i  *  (x] 


' 


-  -V 


!*  aU®  v3  li 


m -l  $&..  .O'--..  ■  'i  {S  ,-rf  '  -  (r;.a'd 


-v' 
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(1?)  b(a,n)  =  b(l,n-a+l)  +  [w(a)  -  w(l)]  b(n-a) 


Using  the  fact  that  b(n,n)  =  g(n)  =  1  ,  we  substitute 
expressions  for  b(a,n)  given  by  (17)  into  b(n)  =  b(l,n)  +  b(2,n)  + 
...  +  b(n,n)  to  obtain 


n-1 


(18)  b(n)  =  1  +  \  b(l,n-a+l)  +  [w(a)  -  w(l)]  b(n-a) 


a=. 


This  can  be  used  to  show 


n-2 


(19) 

b (n)  -  b(n-l)  = 

b(l,n)  + 

\  [w(n-a)  - 

w(n-a-l) ]  b (a )  , 

a^l 

for 

n  >  2  ;  when  n  =  1 

and  2 

we  have 

(20) 

b(i)  .  b(i. 

1)  , 

and  b(2) 

-  b(l)  =  b (1,2)  . 

Equations  (19)  and  (20)  imply  the  following  relationship 
between  the  generating  series l 


00  00  00 


Each  of  the  series  in  (21)  can  be  replaced  with  the  function 
it  represents.  Writing  dB  for  the  partial  derivative  of  B(x,s)  with 


•  (  .  *  ;•  ’  *  i  5 


Ject-ip  vt  <  I  (ais  »  (axa)r  **•  i  "  .a  »/J3 


;  .  -  ;  '  '  ■  ■•  ■■•■■■• 


.  a  o.*  o-u  s'S 


X 


iv-f  h  I  a  nBF.  - 


;  i.  .  ;  c  ;  ,  T)d  *  i 


ft  “*  i$  :.»rf3  OSS^td 


1 

-  n*(a)d 
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respect  to  s  at  s  =  0  ,  the  result  after  collecting  terms  is 


(22) 


SB  =  (w(l)  +  (l-x)(l-W(x)))  B(x,l) 


Now  we  eliminate  dB  from  (12).  First,  setting  a  =  1  in 
(15)  gives 


n-1 


(23) 


b(l,n)  =  w ( 1 )  b(n-l)  +  ^  t(v)  b(v,n-l)  , 

vA 


and  substituting  expressions  for  b(v,n-l)  given  by  (17)  into  the  sum 
in  the  right  member  of  (23)  gives 


n-1 


(24) 


b(l,n)  ac  w(l)  b(n-l)  +  V  t(v)  b(l,n-v) 

v=l 


n-2 


n-2 


^  w(v)  t (v)  b(n-v-l)  -  w(l)  ^  t(v)  b(n-v-l)  l 
v=l  v=l 


for  n  =  1  and  2  ,  the  relations  corresponding  to  (24)  are 


(25) 


b(l, l)  =  1  and  b(l,2)  =  w(l)  b(l)  +  t(l)  . 


Equations  (24)  and  (25)  imply  the  following  relationship 
between  the  generating  functions: 


(26) 


c)B  =  x  +  w(l)x  B(x,l)  +  T(x)  dB 

-  w(l)x  T(x)  B(x,l)  +  xQ(x)  B(x,l)  , 


- 


(S,J)  I 


<  v.  i  .uS  .  j  Vi  *>,.  j.  .  0-1. » :• 

■ 
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where  Q(x)  is  the  function  defined  in  (12) „  Taken  together,  the 
relations  in  (22)  and  (26)  imply  (14). 

Now  to  find  the  generating  function  of  (b(n))  as  defined 
in  (6),  we  put 


(27) 


00 


W(x)  SB  T(x)  ss 

n=l 


/  L  n 
(n-g)x 


x(l+x) 

2(l-x)2 


) 


(28) 


Q(x) 


00 


n=l 


x 


(l+6x+x2) 

4(l-x)5 


Substituting  these  functions  into  (14)  gives 


(29) 


B(x,l) 


x(l-x)? 

2  ^ 

l-5x+7x  -hx 


00 

^  b(n)xn 
n=l 


2  5 

Multiplying  equation  (29)  by  l-5x+7x  -4x^  and  equating 
coefficients  in  the  resulting  identity  gives  b(l)  =  1  ,  b(2)  =  2  , 

b(3)  =  6  ,  b(4)  =  19  ,  and 


(30)  b(n+3)  =  5b(n+2)  -  7b(n+l)  +  4b (n)  , 

for  n  =  2,  3,  ♦ „ „  .  Since  the  largest  real  root  of  the  auxiliary 

equation  for  the  difference  equation  in  (30)  lies  between  3*20  and 
3.21  ,  we  conclude  that  (3*20)n  <  b(n)  <  (3»2l)n  for  sufficiently 


large  n  , 


x  bn* 


-  ^  01  tfoK 

i:  -•  *  < 


'  -  -  r»>  < 


00 


(89) 
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The  method  just  used  to  obtain  a  lower  bound  for  n-ominoes 
also  applies  to  animals  with  cells  of  different  shapes.  For  example, 
animals  with  connected  strips  of  hexagons  in  each  row  are  enumerated  by 


(31) 


b(n)  =  V  (<»1+a2)(a2+a3)  ...  (a. 


where  the  sum  extends  over  all  compositions  of  n  „  This  follows 
since  a  strip  of  r  hexagons  can  be  connected  along  the  upper  edge 
of  a  strip  of  s  hexagons  in  r+s  ways* 

To  find  the  generating  function  for  the  sequence  (b(n)) 
defined  by  (31)  we  substitute 

00 

(32)  w(x)  =  T(x)  =  V  nx11  =  — - — p  ,  and 

U  (1-x)2 


(33) 


00 


Q(x) 


■I 

n=l 


n2xn  _  x£l+xl 


(1-x)- 


> 


into  equation  (l4)  to  obtain 


(34) 


B(x,l)  = 


x 


(1-x)- 


2  ^  k 

l-6x+10x  -7x^+x 


00 

V*  b(n)xn 

n-1 


The  relation  in  (34)  implies  b(l)  =  1,  b(2)  =  3  >  b(3)  =  H  9 

b(4)  =  42  ,  and 

(35)  b(n+4)  *  6b(n+3)  -  10b (n+2)  +  7b(n+l)  -  b(n)  , 


10  8*  1*8  '■  5  iJ  lv  ‘.'Si  OB 


n* 


e+i  j  i  antosaxad  & 


>  -■•  LV  X 

631.  .  -  ■«  )  \4  k  v i:i3b 

7  -  (x)T  -  (x)W 


(«) 


{  a 3 do  o  41)  nciifcupa  oin. 


tJ  «  (l)d  adiiqmi  (4$)  nl  floiixlax  edT 


for  n  =  1,  2,  * . *  o  Furthermore,  the  largest  real  root  of  the  auxiliary 
equation  for  the  difference  equation  in  (35)  lies  between  3*87  and 
3.88  ,  so  that  for  sufficiently  large  n  ,  (3„88)n  >  b(n)  >  (3„8t)R  . 

Golomb  [  ]  suggested  the  problem  of  determining  the  number 

of  incongruent  n  celled  animals  with  hexagonal  cells;  these  numbers 
are  1,  1,  3,  7,  22,  and  83  for  n  =  1,  2,  3,  4,  5,  and  6 
respectively,,  They  correspond  to  free  animals  while  the  numbers  defined 
in  (31)  correspond  to  fixed  animals  of  a  certain  type.  Since  a  hexagon 
has  12  symmetries  in  the  plane,  the  numbers  of  fixed  and  free  animals 
with  hexagonal  cells  differ  by  a  factor  of  12  at  most. 

A  rhombus  is  formed  when  two  equilateral  triangles  are  joined 
along  an  edge;  thus,  a  lower  bound  for  the  number  of  animals  with  n 
rhomboidal  cells  can  be  used  to  find  a  lower  bound  for  the  number  of 
animals  with  2n  triangular  cells „  A  connected  strip  of  r  rhombuses 
can  be  joined  above  a  strip  of  s  rhombuses  in  2 (r+s-l)  ways,  since 
the  strip  of  r  rhombuses  has  two  orientations  with  respect  to  a 
reflection  about  its  midsection*  Thus,  a  lower  bound  for  the  number  of 
fixed  animals  with  n  rhomboidal  cells  is 

(36)  b(n)  =  V  (2a1+2a2-2)(2a2+2a3-2)  ...  (23^+23^2)  , 

where  the  sum  extends  over  all  compositions  of  n  *  Now  to  find  the 
generating  function  of  (b(n))  as  defined  by  (38),  we  can  obtain  the 
appropriate  expression  for  W(x)  and  Q(x)  by  multiplying  equations 


,,*i.  v*  .'io  U«i  J.»***i  o<*.  ,***»«  l««l  •  .  a  -O* 

' 


. 


S  NCI.  ••>  aJ  »:  >0  *  lo  «...  .  ..  «>«*«  oS  bno^etroo  V'  t 


i«,.,  .*  si  ht-otx  *  H  **m*  f****9-1  m*" 


•■  -tv  ■ 


„  yi,  ,»r.-  f  -*)S  fit  «  1  f*»*  *  •***  fc9;,i0t  9j  •» 

c  -  ..  «••©.*  tlatiia  c  r  sad  aoeu  * 


*.!]  ball  03  v#o«  .  n  lo  acoiaieoqsrroo  II«  lavo  *bnd3X9  muz  -u  ®*©riw 


«U  ainldo  nno  6W  ,<*)  xd  b«.iv>b  no  ((«)<*}  io  noldonul  gnll.-woos 


(*>P  bn*  (*)«  *o!  Bi.I^oiqq. 


- 


-  26  - 


(27)  and  (28)  by  2  and  4  respectively.  Substituting  these  functions 
into  (l4)  gives 

(37)  B(x,l)  =  - k  =  V  b(m)xn  . 

l-6x+8x  -6x  -x  L-i 

n=l 

Of  course,  (37)  implies  b(l)  =  1  ,  b(2)  =  3  ,  b(3)  =  13  , 

b(4)  =  59  >  and 

(38)  b(n+4)  =  6b(n+3)  -  8b(n+2)  +  6b(n+l)  +  b(n)  , 

for  n  =  1,  2,  ...  .  The  auxiliary  equation  for  (38)  has  its  largest 

real  root  between  4.5^  and  4.55  >  so  for  sufficiently  large  n  , 
(4.55)n  >  b(n)  >  (4.54)n  ;  from  this  we  can  conclude  that  the  number  of 
fixed  animals  with  n  triangular  cells  is  greater  than  (2„l3)n  for 
sufficiently  large  n  . 
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CHAPTER  III 


IMPROVED  LOWER  BOUNDS  FOR  THE  CELL  GROWTH  PROBLEM 


In  the  last  chapter  we  defined  F  ,  for  each 

^  • 

12  l 

composition  (a^,a^, . . . ,a^)  of  n  ,  as  the  set  of  equivalence  classes 

th 

in  S  containing  n-ominoes  with  a.  cells  in  the  i  row  of  the 
n  j 

n-omino,  for  j  =  1,2,  ...,i  .  We  used  the  fact  that  a  strip  of  s 
cells  can  be  joined  along  the  upper  edge  of  a  strip  of  r  cells  in 
r+s-1  ways;  thus,  there  are  exactly  (a ^ +ag-l ) (a^+a^-l ) . . .  elements 

in  the  subset  F*  of  F  1  containing  n-ominoes  with  a 

a,  a*  . .  .a .  a,a_,««a, 

12  l  1  ^  i 

single  strip  of  cells  in  each  row.  This  gives 


(i) 


f*(al>a2> • • •>ai)  =  f*(a1,a2)f*(a2»a5)-*-f*(al.1)a1)  , 


where  f*p,q,...)  denotes  the  number  of  elements  in  F* 

7  ^ '  pq, 


As  the  relation  in  (l)  suggests,  it  is  easy  to  give  a  one-one 


correspondence  between  the  elements  of  F* 

ai2  *  *  *  i 


and  the  (i-l)-tuples 


of  F*  X  F*  X  ...XF*  .  To  do  this,  suppose  x  is  in  standard 

a, a _  a_a  a.  .a  7 

12  2  3  1-1  i 

position  and  is  an  element  of  the  equivalence  class  X  in  F*  a  a  * 

1  2  *  °  *  i 

The  cells  of  x  in  the  j*”  and  (j+l)St 


rows  above  the  x-axis 


in  a .fMAHo 
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comprise  a  representative  element  y.  of  an  equivalence  class  Y.  in 

J  J 

F*  a  f  for  j  =  1,2, .  ,.,i-l  .  In  this  way,  X  corresponds  to  a 

j  J+l 

unique  (i-l)-tuple  . , , ,Y^  .  On  the  other  hand,  starting 


with  an  element  (Y. ,Yq, . . . ,Y.  ,)  of  F*  X  F*  X  . . . X  F* 

J.  tz.  i  -  j.  a,  a  a_a^  a. 


we 


1  2 


2  3 


i-1  i 


construct  x  as  follows:  Suppose  y^  6  Y^  is  in  standard  position 

and  then  choose  y.  e  Y.  so  that  the  a.  cells  in  the  bottom  row  of 

J  J  J 

Yj  are  the  same  as  the  a^  cells  in  the  top  row  of  y^,  ^  ,  for 

i-1 

j  =  2,3,...,i-l  .  Now  Yj  =  x  is  a  representative  element  of  X 

These  obvious  constructions  establish  a  one-one  correspondence  between 
the  elements  of  the  two  sets  under  discussion. 


Now  we  are  going  to  show  that  in  fact  there  is  a  one-one 


correspondence  between  all  the  elements  of  F  xF  x  ...XF 

ala2  a2a3  ai-iai 


and  a  certain  subset  of  F 


a,  a  . . .a . 

12  l 


Theorem  i;  f  (a^a^  . . .  ,ai)  >  f  (a^a^f  (a^a^) . .  .f  (a^^) 


Proof ;  We  call  an  rxl  rectangle  located  in  a  column  of  the  square 

lattice  and  a  non-empty  subset  of  the  r  cells  contained  in  this 

rectangle  an  r -component .  For  example,  an  n-omino  in  an  equivalence 

class  in  F  is  a  sequence  of  r-components  contained  in  an 

12  r 
r  x  j  rectangle . 

Now  suppose  an  rxj  rectangle  contains  a  sequence  R  of 
r-components  such  that  the  top  row  of  the  rectangle  contains  exactly 
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c  cells  belonging  to  the  components  of  R  ;  also,  suppose  an  s  X  k 
rectangle  contains  a  sequence  S  of  s -components  such  that  the  bottom 
row  of  the  rectangle  contains  exactly  c  cells  belonging  to  the 
components  of  S  .  Clearly,  we  can  translate  the  components  to  the 
right  in  either  or  both  of  the  rectangles  leaving  gaps  between  the 
components  so  that  the  c  cells  in  the  bottom  row  of  S  can  be  made 
to  cover  the  c  cells  in  the  top  row  of  R  .  The  sequence  R*S  of 
(r+s-l) -components  which  results  when  the  components  of  R  and  S 
are  joined  in  this  way  is  called  the  sum  of  R  and  S  . 


X...  X  F 


Let  (Y, ,Y0, . . ,,Y.  ,)  be  a  given  element  of  F  X  F 

jl  d  1-1  a±2  a2a3 

and  suppose  y  is  a  representative  element  of  Y.  , 


a .  .a. 
l-l  i 


j  - r - - j 

for  j  =  l,2,...,i-l  .  We  are  going  to  construct  a  representative 

element  y  of  an  equivalence  class  Y  =  Y(Y. ,Y0, . . .,Y.  « )  of  F 

l  d  ai2  * • #aj 

th  s  t 

such  that  the  sequence  of  2-components  in  the  j  and  (j+l)  rows 
of  y  is  the  same  as  the  sequence  of  2-components  of  y^  ,  for 
J  =  l,2,*o#,l— 1  q 


Consider  the  sequence  of  i-components  given  by 
((...((y1*y2)*y3)*...)*yi_1)  =  yj*y2*. . .*yi<Bl  ;  if  this  sequence  is  not 
an  n-omino,  2-components  containing  disconnected  cells  can  be  translated 
to  the  left  and  joined  to  cells  they  were  joined  to  formerly  in 
y  ,y  , ...,y  ,  it  being  understood  that  overlapping  2-components  must 

Jl  ^  j 

be  translated  simultaneously.  The  n-omino  formed  in  this  way  has  the 
desired  properties;  since,  different  sequences  of 
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Fo  XF  X  ,,,XF  give  rise  to  different  elements  Y  of 

ala2  a2a5  al-lai 

F  ,  the  theorem  is  proved.  An  example  of  the  construction  is 

12  i 

given  in  Figure  1  . 


y4 


yx  yi*y2  (y1*y2)*y,  ((y1*y2)*y5)*y 


Figure  1. 

We  note  that  when  y^y^,**  •  •  ^  (described  in  the  proof  of 

Theorem  l)  is  disconnected,  there  may  be  many  ways  in  which  the  cells 
can  be  translated  to  the  right  or  the  left  to  eventually  form  an  n-omino; 
thus,  n-ominoes  in  different  equivalence  classes  in  F  may 

aiV"ai 

have  the  same  sequence  of  2-components  in  their  jta  and  (j+l)St  rows, 
for  j  =  l,2,...,i-l  .  An  example  of  this  is  given  in  Figure  2  . 
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Both  A  and  B  correspond  to  )  A  corresponds  to  y^  *  y^  . 

Figure  2. 

Since  f(m,  n)  is  generally  larger  than  (m+n-l)  ,  Theorem  1 
leads  to  a  lower  bound  for  s(n)  which  is  larger  than  the  one  given  by 
the  relation  in  (2.6)  .  Thus,  if 


(2)  b(n)  =  ^  f(a1,a2)f(a2,a5)...f(a._1,a1)  , 

where  the  sura  extends  over  all  compositions  (a^,a2, . . . ,a^)  of  n  , 
then  s (n)  >  b(n)  ,  for  n  =  1,2,...  .  Once  we  have  the  generating 
function  of  (f(m,n))  ,  the  methods  given  in  Chapter  I  can  be  applied 
to  the  sum  in  (2)  . 


Theorem  2! 


H(x,y)  =  - 1~^  2 

1-x-y+x  y 


00 


I 

m=o 


00 

f(m,n)xmyn 

n=o 
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Proof: 


First,  we  observe  that 


(4) 


f(n,0)  =  f (0,n)  =  1  ,  n  =  0,1,..., 


and 


(5) 


f(n,l)  =  f(l,n)  =  n  ,  n  =  1,2,.. 


o  t 


Next,  we  show  that 


(6) 


f(m+2,n+2)  =  f(m+2,n+l)  +  f(m+l,n+2)  -  f(m,n)  , 


for  m, n  =  0,1,...  .  There  is  a  way  to  construct  a  representative 

element  of  each  equivalence  class  in  Fm+2  n+2  adding  a  cell  to 

representative  elements  in  each  of  the  equivalence  classes  in 

F  «  .  U  F  ,  .  so  that  exactly  f(m,n)  elements  are  duplicated 

m+2,n+l  m+1, n+2  '  7  ' 

in  the  process.  To  do  this  suppose  x  is  in  standard  position  and  is 

an  element  of  an  equivalence  class  in  F^^  n+i  (or  *n  ^m+l  n+2^  * 

and  join  a  new  cell  to  x  in  the  second  row  (or  first  row)  so  that 

the  new  cell  is  a  maximum  distance  to  the  right  of  the  origin.  It  is 

easy  to  see  that  the  (n+m+4)~ominoes  obtained  in  this  way  are  distinct 

except  those  which  can  be  obtained  by  connecting  a  2x2  block  to  the 

right  end  of  elements  in  standard  position  selected  from  the  equivalence 

classes  of  F  :  these  (m+n+4)-ominoes  will  be  constructed  once  in 

connection  with  an  element  of  F  ,  0  and  again  in  connection  with 

m+i , n+d 


an  element  of  F 


m+2,n+l 


In  order  to  verify  that  H(x,y)  generates  {f(m,n)jm,n  =0,1,...  ) 
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we  check  that  H(x,y)  and  6H(x,y)/8y  respectively  generate  (f(m,0)} 
and  {f(m,  l)}  at  y  =  0;  similarly,  we  see  that  H(x,y)  and 
dH(x,y)/(3x  generate  {f(0,n)}  and  (f(l,n)}  at  x  =  0  .  Multiplying 
relation  (3)  by  (l-x-y+x^y^)  and  equating  coefficients  of  xmyn  ,  we 
see  that  (3)  implies  (6)  .  This  completes  the  proof  of  Theorem  2. 

Using  (4),  (5),  and  (6)  we  can  compute  a  table  of  values  of 
f(m,n)  \  since  f(m,n)  =  f(n,m)  ,  it  is  only  necessary  to  list  f(m,n) 
with  m  >  n  .  The  University  of  Alberta  Computing  Department  found 


for 

m+n  <  30  ,  and  computed 

b(a,n) 

,  a  =  1,2, . . . ,n  ,  and 

n  = 

1,2,..., 30  y  by  means  of 

(1.5); 

these  numbers  are  as  follows 

n 

b(n) 

n 

b(n) 

1 

1 

16 

82830896 

2 

2 

17 

305299746 

3 

6 

18 

1126742091 

4 

19 

19 

4162932345 

5 

63 

20 

15395001510 

6 

216 

21 

56977972065 

7 

756 

22 

211024241503 

8 

2681 

23 

782014829839 

9 

96OO 

24 

2899482522068 

10 

34626 

25 

10755229416160 

11 

125582 

26 

39910544174211 

12 

457^25 

27 

148150485556469 

13 

1671854 

28 

550107559124955 

14 

6127385 

29 

204317^182977793 

15 

22507654 

30 

7590365195768629 

V  * 

.  >r't  to  ‘ 1  30‘u1  <i  Iq  i  f  -i: 

i  •  &  o  a  .03  i  .;j  t  >)  hr  '•  v 

...  o3  xMBioo  n  yliio  «Jt  J 

■'  •  Oil-  -/  '  )  t  »  •  Tj  «;■  .  •  ■'  1  1  f-)  - 


.v  . 
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Theorem  2  gives 

00  00 

(7)  F(x,y)  =  ^  ^  f(m,n)xmyn  =  H(x,y)  -  ; 

m=l  n=l 

furthermore,  we  know  from  the  discussion  in  Chapter  1  that  B(x,y),  the 
generating  function  of  (b(a,n)}  defined  by  (2),  satisfies  the  integral 
equation 

(8)  B(x,y)  =  17^  +  2^;  f  F(x y,  ±)  B(x,s)  f  . 

J  c 

Although  the  kernal  of  the  integral  in  (8)  is  a  rational 
function,  the  solution  function  seems  to  be  a  much  more  complicated 
function,  but  we  have  not  been  able  to  find  B(x,y)  in  closed  form  so 
this  can  only  be  surmised  by  finding  the  functions  B^,Bg,...  described 
in  (1.13),  (1.14),  and  (1.15)  .  The  following  general  technique  can  be 
used  to  estimate  b(n)  from  below.  First,  we  select  a  function 
F*(x,y)  which  generates  (f*(m,n)}  such  that  0  <  f*(m,n)  <  f(m,n)  , 
m,n  =  1,2,.,.,  and  define  a  number  b*(n)  by  substituting  f*  for 
f  in  (2)  ;  of  course,  b*(n)  <  b(n)  .  Furthermore,  we  choose 
F*(x,y)  so  that  it  is  possible  to  obtain  the  solution  function  (in 
closed  form)  of  the  integral  equation 

(9)  B*(x,y)  =  +  2^1  B*(X>8)  T’  ’ 


where  B*(x,y)  generates  (b*(a,n))  . 


-  '  £r 
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/ 
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In  Chapter  I  we  used 

f  f(m,n)  ,  if  m,n  <  j 

(10)  f*(m,n)  =  J 

0  ,  otherwise 

to  define  B^(x,y)  =  B*(x,y)  ,  and  noted  that  B*^(x,l)  is  a  rational 

function  with  its  demominator  equal  to  A^(x)  defined  in  (1.25)  and 

(l,26)  ,  The  polynomials  (-1  )^+^x^  ^+1  ^/2A,  (— )  ,  for  j  =  1,2,..., 6  , 

J  x 

are  as  follows: 

(11)  x-1 

(12)  x3-x2-3x-l 

(13)  x^-x^-3xi+-10x^+2x+l 

(14)  x^0-x^-3x^-10x^-27x^+13x^+l8xi|’+29x^+9x2-5x-2 

(15)  x15-xl4-3x13-10xi2.27x11-70x10+76x9+154x8+242xT+205x6 

+14x3-20x^+120x3-9x2“x+1- 

(16)  x21-x20-3x19-10xl8-27xl7-70xl6-l83x15+35Txli++730x13 

ip  ii  in  Q  8  7 

+1413X  +I743X  +972X  -222x^-665x  -lOOlx1 

-409x8-96x3+378x1++129x3-5x2+x-3 


8 
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From  the  discussion  on  page  9  of  Chapter  I,  we  can  conclude 

that  if  0  is  the  largest  real  root  of  A.(— )  =  0  ,  then 

J  j  x 

b(n)  >  (0.-€)n  ,  for  every  9.  >  €  >  0  ,  and  for  all  sufficiently 
J  J 

large  n  .  Using  (ll)-(l6)  we  find  that  0^  =  1  ,  0^  >  2.41  , 

9^  >  3.00  ,  0^  >  3.34  ,  0^  >  3.50  ,  and  0 6  >  3.59  ;  thus, 


(IT) 


s(n)  >  b(n)  >  (5.59)"  , 


for  all  sufficiently  large  n  . 


There  is  still  another  means  of  approximating  B(x,y)  which 
has  some  advantages  over  the  method  just  used,  but  this  method  depends 
strongly  on  the  form  of  F(x,y)  .  We  can  write 


(18) 


F(x,y)  = 


00 


£=] 


k  V'  fk^x^yk 
Fu(x)y  = 


k= 


bi  (1-*) 


k+1 


where  f^(x)  *-s  a  polynomial  of  degree  M  =  M(k)  .  This  implies  that 

for  fixed  k  and  all  m  >  M  , 


(19) 


k+l 

I 


r=o 


(-Dr  (k!1)  f(m-r,k)  =  0  J 


furthermore,  for  m  >  max  (M(l), . . .  ,M(k)  )  each  of  the  sequences 
(f (m, l) ),  ...,  {f (m,k) )  satisfies  the  difference  equation  (19)  . 


For  a  fixed  j  we  define 


«  0  <  >  <  B  yidvs  iol  *  (»-,  d)  <  a)d 
%  X  *  ^  iftdi  bnil  w  (dl)-(ix)  galeU  .  o  &g*al 


,0'  :  6  ,  Jj.f  <  ,e  ,  oo.f  <  ,• 


(  ?.£  <  (  t)d  <  (o)n 


.  a  ©grrX  yXiaal.Xilua  XIa  loi 


r  (  it  >.  o  q  •  eues  i  ;•  Jons  '  lie  ai  sserfT 

9  jb  »c  h  i  ei::3  Jwd  *  aei;  ja  >oxiiv  '  eril  lavo  89>:  u3navbe  actoe  aari 


■ 


**(*> 


1  i  a  -  ab  lo  x&iraoaylc  a  el  (x)  rl 


- 

*iol  t«ionrr9d3iul 
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e^llsb  .-tv  L9xll  b  to’B. 


-  37  - 


00 


00 


(20) 


where 


F*(x,y)  =  V  F^(x)yk  =  ^  f*(m,n)xmyn 
k=1  m,  n=l 


(21) 


F£(x)  »  ^ 


Fk(x)  ,  for  k  <  M(j) 

j 


V(.l)1+1  (J)  FJ-t(x),  for  k>M(j) 
Cl 


Combining  (20)  and  (21)  we  have 


M  v-1 


I  {  X  (-1)1  Fv-i(x)  } 
(22)  F*(x,y)  =  *=i - £2 - 


(1-y) 


It  follows  from  these  definitions  that 

00  00 

(23)  F(x,y)  -  F*(x,y)  =  V  {f (m,n)-f*(m,n) }  xmyn  , 

m=M+l  n=j+l 

so  that  evidently  F*  tends  to  F  (and  hence  B*  tends  to  B  )  as 
j  tends  to  infinity.  It  also  turns  out  that  0  <  f*(m,n)  <  f(m,n)  , 
for  m,n  =  1,2,..,,  and  j  =  1,2,3  •  This  can  be  checked  by  writing 
F-F*  as  a  rational  function  to  find  a  difference  equation  satisfied 
by  the  numbers  (f (m,n)-f*(m,n) ]  \  the  fact  that  these  numbers  are 
non-negative  can  then  be  established  by  an  induction  argument.  It  may 
be  that  0  <  f*(m,n)  <  f(m,n),  for  j  >  3  ,  however,  we  do  not  need 


this  condition. 


YC  - 


00 


V(*)^ 


s  V;  :.ti  . 


»pc  b  b  1:  n  it  Ii  no/Jai  a  »»  *1- i 
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If  we  put  j  =  3  in  (22)  and  use  the  resulting  function 
F*(x,y)  as  the  kernal  in 

(24)  B*(X,y)  =  +  ITT  /  F*(xy’  7)  B*(X>s)  T  ’ 

J  c 

the  integral  equation  becomes 


(25) 


B*(x,y)  “  1^7  +  {3F1(xy)-3F2(xy)+F5(Xy))B*(x,l) 

+  {-2F1(xy)+3F£(xy)-F3(xy))  dB*(x,l) 
+  {F1(xy)-2F2(xy)+F5(xy)}  d.Bfof..1)  , 


where  d  B(x, l)  denotes  d  B(x,s)/ds  evaluated  at  s  =  1  .  Letting 

Q^(x,y)  ,  Q2(x,y)  ,  and  Q^(x,y)  denote  the  coefficients  of  B*(x,l), 

o 

dB*(x, l)  ,  and  d  B*(x, l)/2  respectively,  we  take  partial  derivatives 
with  respect  to  y  at  y  =  1  in  (25)  to  obtain  a  system  of  equations 
equivalent  to  the  following: 


2 

(26)  {Q1(x,l)-l}B*(x,l)^2(x,l)aB*(x,l)+Q3(x,l) 


(27) 


-X 


(l-x)‘ 


2 

SQ1(x,l)B*(x,l)+(dQ2(x,l)-l)aB*(x,l)-^Q5(x,l)  5 


(28) 


2  8  Q.  (x,l)  3  Q£(x,1)  r8  Q,(x,l)  'i  ^2g#/ 

1  b*(x,i)+— ^ - 8B*(x,i)+-|  — 2_ - 1| 


(1-x)- 
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Applying  Cramer's  rule  to  this  system  to  solve  for  B*(x,l), 
we  obtain  a  rational  function  having  its  denominator  equal  to 
A(x)(l-x)12  where 


F1  -  1 

F2  -  1 

F  -  1 

3 

Sf.  -  1 

dF0  -  2 

dFx  -  3 

(29) 

A(X)  = 

1 

2 

3 

a2Fi 

1 

•ri 

°io  1 

2 

2  1 

2 

i  dLF,(xy) 
and  d  F.  =  - ^7 - 

j  Sy1 


.  Substituting  the  values  of 


y=i 


1  P 

given  by  (7),  we  expand  (29)  to  find  that  A(x)(l-x)  is  equal  to 


(30)  3x^-15x*2-f^3x^-1^3xl^+288x9-662x^+1101x^-126lx^+1032x^ 

-6l2xi++260x^-75x2+13x-l 

12 

The  largest  real  root  of  A(l/x) (l-l/x)  =0  is  between 

3.60  and  3.61;  hence, 

(31)  s(n)  >  b(n)  >  b*(n)  >  ( 3 . 60 ) n  , 


for  all  sufficiently  large  n  . 

If  regular  hexagons  of  unit  area  are  used  as  cells  of  an 

animal,  the  number  of  ’’two  row  animals"  (analogous  to  the  elements  of 

F  )  of  this  type  with  m  cells  in  the  first  row  and  n  cells  in  the 
mn  '  J  r 


■'V 


kA 


(oO 

j  ^rfX 

.  io» '  :i  ;  f  f  kftfi  0  .$. 
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second  is  exactly  (m+n)  .  By  modifying  the  proof  of  Theorem  I  slightly, 
we  can  show  that  there  are  more  than 


(32) 


b(n) 


a  +a 
12 


a. 


a  +a 

2  3 


ai-i+ai 


’’fixed"  animals  with  n  hexagonal  cells,  where  the  sum  in  (32)  extends 
over  all  compositions  (a^,a^, . . . ,a^)  of  n  *  It  follows  from  (5*3) 
and  (5.10)  that  b(n)  >  p(n+2)  =  (^n)  .  Thus,  for  every  0  <  c  <  4  , 


(33)  b(n)  >  (It.  -  «)n  , 

for  all  sufficiently  large  n  . 


*  i  0^  (Oi*  )  -;n'• 
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CHAPTER  IV 


READ'S  THEOREM 


We  mentioned  in  Chapter  II  that  Read  pL7  ]  enumerated  the 
number  of  equivalence  classes  in  containing  n-ominoes  with  cells 

in  exactly  k  rows  of  the  plane.  He  characterized  such  an  n-omino  as 
a  sequence  of  k-components  which  satisifies  certain  restrictions:  Some 
k-components  are  permitted  to  begin  a  sequence >  each  k-component  i 
may  be  succeeded  in  the  sequence  only  by  k-components  selected  from  a 
set  X^,  and  only  certain  k-components  may  terminate  a  sequence.  In 
what  follows,  we  reformulate  the  problem  which  Read  posed  and  solved, 
and  give  an  alternative  solution  based  on  the  methods  outlined  in 
Chapter  I, 

Let  X  =  {1,  2,  ...,  k}  be  a  given  finite  set  of  objects  and 

suppose  Y  and  Z  are  non-empty  subsets  of  X  .  Furthermore,  suppose 

that  each  object  i  in  X  has  assigned  to  it  a  subset  X.^  of  X  and 

an  integer  w(i)  .  We  let  t(q,n)  denote  the  number  of  q-tuples 

(a1,a2,...,a  ),  a±  €  X  ,  i  =  1,  2,  ...,  q  ,  such  that 

(i)  w(ax)  +  ...  +  w(aq)  =  n  ,  (ii)  a^  €  Y  ,  (iii)  aq  e  Z  , 

and  (iv)  a.  €  X  ,  i  =  2,  3,  q  .  Now 

1  2  .  .. 


-N 


,1  I  •!  •  l  ■'  j*  *  j  I  •  O'  1  ins®  •'« 


■jM .  ...  1  V)  3=.:»-psa  » 


- 


t  co  bt,a«  noil  o  i  >v  t3*c toil*  n&  v  g  Dn6 
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(1)  Tq(z)  =  ^  t(q,n)zn  ,  and 

n 

^  t(q,n)sqzn 

q  n 

th 

are  called  the  q _  counting  polynomial  and  the  sequence  counting  series 

respectively.  Read  proved  the  following  theorem  which  gives  a  means  of 
calculating  T^(z)  and  T(s,z)  . 

Theorem  (R.  C.  Read): 


(3) 


T  (z)  =  P[V(z)]q'1S  ,  and 


00 


-1. 


T(s,z)  =  sP[I-s^(z)]“  S  , 


where  I  is  the  k  X  k  identity  matrix,  S  is  the  column  vector 
{s^Sg,  ...,sk]  ,  P  is  the  row  vector  [Pl>P2>  •  •  *>Pk)  >  and  ^(z) 
is  a  k  X  k  matrix  [f(m, n) ]  defined  by 


(5) 


i  a  4 

11  l 


w(n) 


if  n  e  Y 


o  ,  otherwise 


1  ,  if  n  e  Z 


n 


(6) 


0 


t 


otherwise 


and 


s  a{n,p)j  ^  V 


( 


3  a  ilc  >  .*  ?  *r  ,x  -a  a*  y:*Jbi9bl  X  >!  oril  il  I  sxsrfv.' 


■ 


,rt  'Ws 

o 

, 

ty-  ■■ 
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(7) 


if  n  e  X 


m 


otherwise 


We  will  not  present  Read's  proof  of  his  theorem,  but  proceed 
at  once  in  deriving  a  second  solution  of  his  problem.  First,  using  the 
definitions  of  f(m, n)  and  p^  given  in  (6)  and  (7)  we  define 


(8) 


F(x,y:z) 


f  (m,  n) 


m  n 

x  y  > 


and 


(9) 


f 


and  then  put  sF(x,y:z)  and  sG(x)  in  place  of  F(x,y)  and  G(x) 
in  (1.7)  to  define  b(n)  as  in  (l.l)  .  Now  consider  the  generating 
function  of  the  sequence  {b(a,n)}  defined  in  this  way: 


(10) 


00 

B(x,y:s,z)  = 

n=l 


If  (a^,a  ,  ...,a  )  is  a  composition  of  n  which  satisfies  the  four 

,,rj  ;•  i’ 

properties  used  to  define  t(q,n)  ,  then  it  contributes 

a. 


w(a2)+w(a  )+...+w(a  ) 
sqz  ^  q 


1  n 

to  the  coefficient  of  y  x  in  (10); 


otherwise,  the  contribution  is  0  .  We  want  the  sum  of  the  generating 
functions  of  {zW^b(a,n) :  n  =  1,2,...  ]  for  a  e  Y  ;  thus, 


(11) 


T(s,z)  =  V  ^7  ^  (B(l>y-S>z)} 

n»l  *  ^ 


y=o 


•  .» 


) 

. 


s-v  8 


{  . . n  :(ctr )  ^ 
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We  use  the  relations  in  (1.24),  (1.25),  and  (1.26)  to  find 

the  partial  derivatives  in  (ll).  Let  <X>(s,z)  =  [sf(m,n)],  let  I 

denote  the  identity  matrix  of  order  k  ,  let  A  =  ®(s,z)  -  I  ,  and 

th 

let  A^  denote  the  matrix  obtained  by  replacing  the  n  column  of 

A  with  the  column  vector  {-sp^,  -sp^,  . ..,  -sp^}  •  This  gives 

k 

(12)  T(s,z)  =  ^  s^  det  A^/det  A 

The  following  example  illustrates  the  use  of  Read’s  Theorem 
and  relation  (12).  Letting  1,  2,  and  3  denote  the  2-components  with 
one  cell  in  the  top  row,  with  one  cell  in  the  bottom  row,  and  with 
two  cells  respectively,  we  put  w(l)  =  w(2)  =  1  ,  w(3)  =  2  ,  =  {1,3}  , 

X2  =  {2,3}  ,  X^  =  Y  =  Z  =  {1,2,3}  t  now  t (q, n)  is  the  number  of 
n-ominoes  which  may  be  constructed  as  a  sequence  of  q  2-components. 

Read  used  his  theorem  to  show  that  in  this  case 


sz-1 

0 

2  “1 

sz 

[-1 

z 

(13) 

T(s,z)  =  [ s , s , s  ] 
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Using  relation  (12)  we  find 
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sz-l 

0 

c. 

SZ 

(14) 

det  A  = 

0 

sz-l 
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sz 

sz 

sz 

2  . 

sz  -1 

=  (sz-l) (l-sz-sz^-s^z^)  , 


(15) 


(16) 


(17) 


det  A^  = 


det  A2  = 


det  A_  = 
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-s 

0 

sz-l 

-s 

sz 

sz 

-s 

=  s(sz-l)  , 


=  s(sz-l) 


=  s  (sz-l) (sz+l) 


so  that  combining  these  we  obtain 


(18) 


2 

T(s,z)  «  (z  det  A^  +  z  det  Ag  +  z  det  A^]/det 


o  2  2  3 

2sz  +  sz  +  s  z 

.  2  2  3* 
1  -  sz  -  sz  -  s  z 


The  method  outlined  in  Chapter  I  can  also  be  used  to  find 
sequences  (t*(q,n):  q,n  =  1,2,...  }  and  (t^(q,n):  q,r  =  1,2,...}  , 
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such  that 


(19) 

t*(q,n) 

<  t(q,n)  <  t*(q,n) 

y 

for  q, 

n  =  1,2, ... 

To  do  this  we  let 

denote  the  subset  of  X 

containing  all  objects 

with  weight  i  ,  and  define 

(20) 

f*(m,  n) 

=  |x  n  w  |  , 

j  €  w  3 

J  m 

and 

(21) 

f*(m,n) 

=  min  |X  f|  W  1  , 

j  6  W  J  ” 

m 

where 

| S |  denotes  the 

i  number  of  elements  in 

set  S  .  Also,  we  define 

(22) 

g*(n)  = 

J  1,  if  wn  n  Z  4  <} 

^  0,  otherwise 

,  and 

(25) 

g*(n)  = 

r  1,  if  w  c  z 
{ 

^  0,  otherwise 

and  let 

(24) 

b*(a,n)  =  ^ 

s1  f*(ax,a2)  f*(a2,a^) 

...  g*(ai)  > 

(25) 

b*(a,n)  =  ^ 

...  (a  j^^^y a^ )  g^.(aj^)  y 

where  the  sums  in  (24)  and  (25)  extend  over  all  compositions  (a^a  ,  . .  .#3^ 
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of  n  with  =  a  .  The  generating  functions 


(26) 


B*(x,yts)  =  b*(a,n)yaxn  ,  and 


a,n 


(27) 


are  defined  by 


B*(x,y:s)  =  ^  b^(a,n)yaxn  , 


a.n 


(28) 


sF*(x,y)  =  ^  f*(m,n)sxmyn  ,  sG*(x)  =  ^ g*(n)sxn  ,  and 


m,n 


n 


(29) 


sF*(x,y)  =  ^  f*(m,n)sxmyn  ,  sG*(x)  =  g*(n) 


sx 


n 


m,  n  n 

respectively,  just  as  B(x,y)  is  defined  by  F(x,y)  and  G(x)  in 
(l.l).  Letting  s(a)  denote  the  number  of  objects  in  Y  with  weight 
a  ,  we  define 


(50) 


T*(s,z)  =  ^ 


s(ay 

a  l 


rs  (»*(*, y;s))0 

oy 


and 


(31) 


(B,(z,y:s))y=0 


The  sequences  (t*(q,n))  and  {t*(q,n))  generated  by  (30) 
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and  (31)  respectively,  satisfy  (19);  this  follows  from  (24)  and  (25), 
and  the  definitions  of  f*,  g*,  f#,  and  .  Suppose  (a^,a2, . . . ,a^) 
is  a  composition  of  n  and  consider  its  contribution  to  the  sums  in 
(24)  and  (25):  There  are  at  most  s(a1)f*(a1,a2). . ,f*(aq  1>aq)§*(aq) 
and  at  least  s  .  •  •f*(aq_1,a^)g^(a^)  sequences 

of  objects  (O^Og,... ,0  )  such  that  w(0.)  =  a.  ,  for  i  =  1,2,  ...,q  „ 
The  coefficient  of  s^z"  in  (JO)  and  (31)  is 


(52) 


and 


t*(q,n)  =  Vs^)  f*(a2,a^)  ...  f*(a^_1,a(])  g*(a(])  , 


(33) 


t*(q,n)  =  2^8(8^  f*(a1>a2)  f,(a2,a3)  ...  f*(aq.i>aq)  g*(aq)  , 


respectively,  where  the  sums  in  (32)  and  (33)  extend  over  all  compositions 
of  n  into  exactly  q  parts. 

The  sums  in  (30)  and  (31)  are  finite,  and  the  partial  derivatives 
involved  can  be  found  in  terms  of  determinants  as  shown  in  (1.24)  -  (1.26). 
From  this  we  can  conclude  that  T*(s,z)  and  T^(s,z)  are  rational 
functions  just  as  T(s,z)  is.  Furthermore,  if  f*(m, n)  =  f^(m, n) 
and  g*(n)  =  g^(n)  ,  for  m,n  =  1,2,...,  then  T*(s,z)  =  T^(s,z)  =  T(s,z). 
This  last  fact  can  sometimes  be  used  to  greatly  reduce  the  work  involved 
in  calculating  T(s,z)  by  means  of  Read’s  Theorem  or  the  formula  given 
in  (12). 
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The  example  treated  in  (13)  -  (18)  is  a  case  for  which 
T*(s,z)  =  T#(s,z)  =  T(s,z):  This  follows  since  all  of  the  objects  are 
weighted  1  or  2  ,  and  f*(l,l)  =  f*(l,l)  =  1  ,  f*(l,2)  =  f*(l,2)  =  1  , 
f*(2,l)  =  f^(2,l)  =  2  ,  f*(2,2)  =  f^(2,2)  =  1  ,  and  any  of  the  objects 
may  begin  or  end  a  sequence.  Calculating  T(s,z)  by  means  of  (30)  and 
(1.24)  -  (1.26)  we  obtain: 


(34) 


T(s,z) 


-sz 

sz 

sz-1 

-sz 
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2  , 

+ 

0  2 
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-sz 

sz  -1 

2sz 
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sz-1  sz 


o  2  2  1 

2sz  sz  -1 


O  2  2  3 

2sz  +  sz  +  s  z 

.  2  2  3 
1  -  sz  -  sz  -  s  z 


Before  we  leave  this  example  we  recall  a  result  implied  by 
Theorem  2  of  Chapter  III;  If  H(x,y)  =  (l+xy )/(l-x-y+x  y  )  generates 
{f(m, n))  ,  where  f(m, n)  denotes  the  number  of  (m+n) -ominoes  with 
m  cells  in  the  first  row  and  n  cells  in  the  second  row,  then 
H(x,x)-1  generates  the  number  of  two  row  v-ominoes,  for  v  «  1,2,... 
Since  T(l,x)  also  generates  the  number  of  two  row  v-ominoes,  for 
v  1,2,...,  we  must  have  T(l,x)  =  H(x,x)-1  ;  the  fact  that  this  is 
indeed  the  case  provides  a  check  for  our  work. 

If  we  want  to  obtain  upper  and  lower  bounds  for  the  numbers 
t(n)  we  use  the  fact  that 


(35) 


t*(n)  <  t(“)  <  t*(n) 
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which  follows  after  summing  on  q  in  (19).  It  has  already  been  observed 
that  the  sequences  {t^(n)}  ,  {t(n)j  ,  and  (t*(n))  are  generated  by 
the  rational  functions  1^(1, z)  ,  T(l,z)  ,  and  T*(l,z)  respectively. 
If  the  sequences  are  increasing. 


(36) 

lim  (t*(n))Vn 

=  (X  , 

lim  (tCn))1/11  =  P 

,  and 

n  -+00 

n  ->  00 

lim 

(t*(n))Vn 

=  7 

n  -*  00 

all  exist 

and  are  equal  to 

the  largest  real  roots  of 

the  auxilliary 

equations 

associated  with 

,  T(1 

, z )  ,  and 

T*(l,z)  respectively. 

For  example,  the  equation 

associated  with 

T*(l,z) 

is 

(37) 

det  [a  ] 
u  mn 

=  0 

y 

where 

r  f*(m. 

m)-zm  , 

if  m  as 

n 

(38) 

a  = 
mn 

a.  f*(m, 

n)  » 

if  m  ^ 

• 

y 

n 

the  equation  associated  with  T^(l,z)  is  defined  by  writing  the 
asterisks  in  (38)  as  subscripts  instead  of  superscripts. 

If  there  are  two  or  more  objects  with  the  same  weight,  the 
determinants  used  to  calculate  T*  or  T*  are  smaller  than  those  used 
to  find  T  .  If  relations  (37)  and  (38)  are  used  to  find  a  and  7  , 
it  follows  from  (35)  that  for  given  €  ,  a  >  e  >  0  , 
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(59)  (a-e)n  <  t*(n)  <  t(n)  <  t*(n)  <  (?+«)"  , 

for  all  sufficiently  large  n  .  Furthermore,  (39)  provides  sharp 
bounds  for  t(n)  if  f*(m,n)  =  f*(m,n)  ,  and  g*(n)  =  g*(n)  ,  for 
m, n  =  1,2,...,  since  then  a  =  (3  =  y  ,  Note  also  that  if  there  is 
no  more  than  one  object  with  a  given  weight  then  f^(m, n)  =  f(m, n) 

=  f*(m,n)  ,  and  g*(n)  =  Pn  =  g*(n)  .  Also,  we  might  expect  to 
obtain  "good"  bounds  for  t(n)  if  there  is  only  "slight"  variation  in 
the  sets  of  numbers  used  to  define  f*(m,n)  or  f^(m,n)  in  (20)  and 
(21). 

The  sequence  (t(n)}  generated  by  T(l,z)  in  (34)  gives 

the  number  of  two  row  n-ominoes  for  n  =  1,2,...  .  In  this  case  a  , 

3  ,  and  7  as  defined  in  (3 6)  are  all  equal  to  the  largest  real  root 

3  2 

of  z  -z  -z-1  =  0  which  is  between  1.84  and  1,85  .  Thus,  we  can 
conclude  that  if  t(n)  denotes  the  number  of  two  row  n-ominoes  then 

(40)  (1.84)"  <  t(n)  <  (1.85)"  , 

for  all  sufficiently  large  n  . 

Read  constructed  the  three  row  n-ominoes  using  the  components 
illustrated  in  Figure  1(a)  as  the  objects  in  his  theorem;  the  label  and 
object  number  of  each  component  is  also  indicated  in  the  figure.  Unlike 
the  two  row  case,  the  first  k  terms  of  a  sequence  of  components 
corresponding  to  a  three  row  n-omino  need  not  form  a  connected  set; 
such  an  n-omino  is  illustrated  in  Figure  l(b). 
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12  3k  5678 


Figure  1(a) 


The  weight  of  a  component  is  the  number  of  cells  it  contains af 
we  put  X  =  (1,2, 3, 4, 5, 6, 7, 8)  ,  Y  =  X  -  (4)  ,  Z  =  X  -  {5}  , 
xx  =  x  -  (5}  ,  Xg  =  x  -  (4,7 )  ,  x5  =  x  -  (4,6)  ,  x4  =  x  -  (5,8)  , 
x5  =  (1,5)  ,  x6  =  (1,2,5, 6)  ,  x?  =  (1, 3,5,7)  ,  xg  =  {1,2,3,83  , 

and  apply  Read's  Theorem  to  obtain  a  generating  function  for  the 
number  of  3X  q  n-ominoes: 
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T(b,z)  =  sP[I  -  sy(z)]  s  , 
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P  =  row  vector  [1,1, 1,1, 0,1, 1,1]  , 


(43) 
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S  =  column  vector  [z  ,z  ,z  ,0,z  ,z,z,z]  , 
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we  would  have 

to 

find 

largest  real  root  of  det  [I  -  ) ]  =  0  \  Instead  of  doing  this,  we  will 


estimate  3  by  means  of  (39)  »  We  have  under  present  assumptions 
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12  1 

(45)  [ f*(m, n) ]  = 

2  3  1 

3  3  1 

,  and  [ f*(m, n) ]  = 

Oil 

3  3  1 

so  that  the  auxiliary  equations  associated  with  T*(l,z)  and  T^(l,z) 

6*54^2  654^2 

are  x  -x  -3x  -2x  -2x  +1=0,  and  x  -x  -x  -2x  -2x  -1=0, 

respectively.  The  largest  real  root  of  the  first  of  these  equations  is 

between  2.57  and  2.58  and  the  largest  real  root  of  the  second  is  between 

2.09  and  2.10  so  that  if  t(n)  denotes  the  number  of  three  row  n-ominoes, 

then 
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(2.0  9)"  <  t  (n)  <  (2.58)" 


for  all  sufficiently  large  n  . 
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We  conclude  this  chapter  by  finding  difference  equations  for 
the  number  of  three  row  n-ominoes,  and  at  the  same  time  illustrate 
another  enumeration  technique  related  to  Read's  Theorem. 

If  X+m+n  =  V  ,  we  again  let  F.  denote  the  set  of 

7  imn 

equivalence  classes  in  containing  v-ominoes  with  i,.m,  and  n 

cells  in  the  first,  second,  and  third  rows  respectively;  also, 

f(i,m,n)  is  the  number  of  elements  in  F .  ,  If  V  denotes  a 

component  illustrated  in  Figure  l(a),  we  let  denote  the  subset 

of  equivalence  classes  in  F^,mn  containing  v-ominoes  whose  sequences 

of  components  end  with  V  ;  also,  v(i,m, n)  is  the  number  of  elements 

in  and  V(x,y,z)  is  the  generating  function  of  the  sequence 

[v(i,m,n):  i,m,n  =  0,1,2,...  )  .  We  write  (V  U  W  U  for 

V„  U  W  U  ...  ,  and  (v+w+. . . ) (i,m, n)  for  the  number  of  elements 
iron  imn 

in  this  set. 

The  component  C*  cannot  end  a  sequence  of  components  which 

form  a  v-omino,  so  C*  is  empty;  thus,  we  have  by  definition 

imn 

(1+7 )  (f*-a-b-b*-c-d-d*-e)  (i,m,n)  =  0  , 

for  i,m, n  =0,1,...  ,  where  F^^^  =  (A  U  B  U  ^mn  .  We  will  see 
presently  that  F^  is  contained  in  F*mn  . 

Also,  we  have  the  following  relations? 

(48)  b(i,m,n)  -  (f*-d*) (i-l,m-l,n)  =  0  ,  i,m  >  1  ,  n  >  0  , 


,  Ri  d  s*  .  i  U  yd  ^qarfo  elrLf  abuioi  oa  sV 


,1  •  ' 


.i  »io3ffT  a'ba&Ji  oi  b*3*Iai  aupiflitoaS  noliMxnuns  lorfaous 
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j  a^u-^ala  io  n  » ’1  «1  (nKnixi  ri: 
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of  ,  oalroo-v  a  oaol 

r  f  ’  *  ' 


!  ' 


••:  w"  '•  •  v.^nsaaiq 
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m,  n  >  1  , 


(49) 

b*(i,m,n) 

-  (f*-d)(^,m-l,n-l)  =  0  , 

(50) 

c(i,m,n)  - 

(a+c) (i-l,m,n-l)  =  0  , 

(51) 

d(^,m, n)  - 

(a+b+c+d ) (^-l,m, n)  =  0  , 

(52) 

d*(i,m,n) 

-  (a+b*+c+d*) (i,m, n-l)  =  0 

(53) 

e(i,m,n)  - 

(a+b+b*-ie)(i:,m-l,n)  =  0  , 

£  >  0 


m  >  0  ,  £,n  >  1  , 


m,  n  >  0  ,  £  >  1 


,  £,m  >  0  ,  n  > 


i,n  >  0  ,  m  > 


> 


1 


1 


To  prove  (48),  note  that  if  B  is  deleted  from  each  v-omino 
in  B^,mn  ,  £,m  >  1  f  n  >  0  ,  then  each  of  the  elements  of  the 
equivalence  classes  in 

(54)  (AUBUB*UCUDUE)... 

<B-l,m-l,n 

is  obtained  exactly  once,  but  this  is  the  set  (F*-B*)  .  .  .  The 

A/"  1  ^  Ttl03  JL  y  n. 

existence  of  this  one-one  correspondence  between  the  elements  of  B  . 

imn 

and  those  of  (f*-D*).  .  1  implies  (48);  the  proofs  of  (49,  -  (53) 

lu~  l ,  m- 1 ,  n 

are  similar. 

If  A  is  deleted  from  each  v-omino  in  A  .  we  obtain  the 

imn 

elements  of  F*  „  ,  ,  and  some  disconnected  sets  of  cells.  The 

i-l,m-l,n-l 

v-ominoes  in  A.  which  give  rise  to  disconnected  sets  of  cells  may 

imn 

be  characterized  as  sequences  of  components  in  which  the  last  k+2 
terms  have  one  of  the  forms  ( . . . ,B,C*, . . . ,C*, A)  ,  ( . . . ,B*,C*, . . . ,C*,A) 

/  St 


*  '  4>’  i*. 


*  C  f 


'  ft  ,rr.  . )  > 
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(. . .,D,C*, .. ,,C*,A)  ,  or  (...,D*,C*, ...,C*,A)  ,  for  k=  1,2,...  ; 

we  also  have  a  sequence  (D, . . . ,D,C*, . . . ,C*,A)  if  Jt,  >  n  ,  (C*,  . . . ,C*, A) 

if  H  =  n  ,  and  (D*, , . . , D*,C*, . . . ,C*, A)  if  i  <  n  .  The  number  of 

elements  in  A„  is  therefore, 

imn  f 


(55) 


a(i,m,n)  = 


l+f*(i-l,m-l, n-l) 


+ 


00 

Z(b+b*+d+d*)  (i-k,m-l,n-k) 
».=2 


Using  the  last  expression  for  a(i,m, n)  ,  we  find  the 
difference  a(^,m,n)  -  a(i-l,m,n-l)  ,  and  obtain  after  transposing  a 
term, 


(56)  a(i,m,n)  =  a(^-l,m,n-l)  +  f*(i-l,m-l,n-l) 

+  (b+b*+d+d*-f*)(i-2,m-l,n-2) 

=  2(^-1,  m,  n-l)  +  f*(j0-l,m-l,n-l) 
-  (a+c+e) (i-2,m-l,n- 2)  , 


for  i!,n  >  2  ,  m  >  1  . 

Each  element  of  B  „  or  D.  can  be  reflected  about  its 

^mn  iimn 

midsection  to  obtain  elements  of  B*^  or  D*m^  respectively,  and 
vice  versa;  from  this  it  follows  that 

(57)  b(^,m,n)  =  b*(n,m,i),  and  d(i,m,n)  =  d*(n,m,i)  . 


The  difference  equations  in  (48)  -  (53) >  (56)  and  (57) 


1 


. 


i  .  ♦  ,  J(  *€  l t bnii  K  Cl  3l  IX 


I  .  (n -a  li)i 


;• 
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provide  a  means  of  calculating  f*(,0,m,  n)  without  using  generating 
series.  The  sets  F*  differ  from  the  sets  F  in  that  v-ominoes  in 
F  have  cells  in  exactly  three  rows,  but  this  is  not  necessarily  the 
case  for  v-ominoes  in  F*  .  These  difference  equations  imply  relations 
between  the  generating  series? 


(58) 


(F*-A-B-B*-C-D-D*-E)(x,y,z)  =  0  , 


(59) 


(B-xyF*+xyD)(x,y,z)  =  xy  , 


(60) 


(B*-yzF*+yzD*)(x,y,z)  =  yz  , 


(61) 


(C-  xzA-xzC  )  (x,  y  ,  z  ) 


xz  , 


(62) 


(D-xA-xB-xC-xD) (x,y, z) 


(63) 


(D*-zA-zB*-zC-zD*) (x,y, z) 


z 


(64) 


(E-yA-yB-yB*-yE)(x,y,z)  =  y  , 


(65) 


(66) 


B(x,y,z)  =  B*(z,y,x)  , 


(67) 


D(x,y,z)  =  D*(z,y,x) 
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When  x  =  y  =  z  ,  this  system  reduces  to  a  system  of  six 
equations  linear  in  the  functions  F*,  A,  B  =  B*,  C,  D  =  D*,  and  E  . 
Applying  Cramer's  rule  to  this  system,  we  are  led  at  once  to  a  determinant 
of  a  matrix  having  the  form  of  (kk)  . 


u  J[o  :  •.  8' '  f;  S  oi  *  'J  i  s  *  K  *  x  t  3<lW 

. 

3  d  o  3  4J  3  ^  a  OJ  *-  L  H*1  •  g  l^IqqA 


-  59  - 


CHAPTER  V 


PLANE  TREES 


In  the  proof  of  Lemma  3  in  Chapter  II  we  described  Eden's 
method  for  assigning  a  sequence  of  binary  digits,  W(X)  ,  to  each 
n-omino,  X  .  If  W(X)  =  (p^, p^, . . . ,  P^n  ,  we  must  have 

=  ^3n-2  =  ^3n-l  ~  ®  ’  since  Cn  cannot  have  any  directed  edges 
going  from  it  to  other  cells  in  X  .  Thus,  the  sequence  (0,  P1,pg, . . 

P^n-4^  =  ^al>CL2f  * '  *  is  SUch  that 


(1) 


33 


1=] 


ai  >  j  , 


for  j  =  1,2, ...,n-l,  since  (0,P^,Pg, . . „,p^, 0,0,0, ... )  corresponds 
to  an  r-omino  (r  >  j)  which  forms  a  part  of  X  .  We  will  show 
presently  that  there  is  a  natural  one-one  correspondence  between  binary 
sequences  (a^,a^, . . . ,a^n  *)  which  satisfy  (l),  and  the  elements  of  a 
class  of  trees  embedded  in  the  Euclidean  plane  having  nodes  with  valency 
1  or  4  .  Besides  exploring  this  idea  further,  we  will  show  that  various 
kinds  of  plane  trees  are  enumerated  by  sums  of  the  type  studied  in 
Chapter  I  . 


Harary,  Prins,  and  Tutte  have  written  joint  and  separate  papers 
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on  the  subject  of  plane  trees  [£],  [7  ],  £>3];  for  completeness  we 
reproduce  their  definitions  here,  A  tree  is  a  connected  graph  which 
contains  no  polygons,  and  a  plane  tree  is  a  tree  embedded  in  the  plane 
with  no  distinct  edges  intersecting,  A  rooted  (plane)  tree  is  a  (plane) 
tree  in  which  one  node  is  distinguished  as  a  root .  A  planted  plane  tree 
is  a  rooted  plane  tree  in  which  the  root  has  valence  1  0  Two  trees 
are  isomorphic  if  there  is  a  one-one  mapping  of  the  nodes  of  one  onto 
the  nodes  of  the  other  which  preserves  adjacency.  Two  plane  trees  are 
map-isotnorphic  if  there  exists  an  orientation  preserving  homeomorphism 
of  the  plane  onto  itself  which  maps  one  onto  the  other*  Similarly,  two 
rooted  (plane)  trees  are  (map~)  isomorphic  if  there  exists  a  (map-) 
isomorphism  between  them  in  which  the  roots  correspond. 


We  let  denote  the  set  of  non-map- isomorphic,  planted, 

plane  trees  with  n  nodes,  and  let  p(n)  denote  the  number  of  elements 

in  P  ,  The  level  of  a  node  in  an  element  of  P  0  is  the  length  of 
n  n+£. 

the  path  from  the  root  to  that  node.  Corresponding  to  each  composition 

(l,  l,ag,a  ,  0 . . )  of  n-j-2  into  an  unrestricted  number  of  positive  parts, 

there  is  a  subset  of  P  ^  containing  elements  which  have  exactly  a. 

n+2  1 

nodes  with  level  i  ,  for  i  =  2,3,...  .  The  nodes  having  the  same 
level  can  be  arranged  in  rows  so  that  no  edges  in  the  tree  cross,  and 
the  nodes  in  each  row  can  be  labelled  consecutively  from  left  to  right 
in  successive  rows  beginning  with  the  root  of  the  tree. 


Evidently,  all  of  the  edges  of  a  tree  in 

til 

consecutive  levels  such  that  each  node  in  the  j 


P  „  occur  between 
n+2 

level  is  connected 
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to  exactly  one  node  in  the  (j-l)St  level,  and  every  node  in  the  (j-l) 

level  is  connected  to  a  (possibly  empty)  set  of  consecutively  labelled 

th 

nodes  in  the  j  level.  If  a.  ,  and  a.  denote  the  number  of  nodes 

J-l  J 

S  t 

in  the  (j-l)  and  j  levels  respectively,  then  the  number  of  ways 


st 


edges  may  be  drawn  between  these  levels  is 


.  /aj-i + 

IS  j 

V  aj 


-i +  aj  •  1 


the  number 


of  ways  of  selecting  a.  objects  from  a  set  of  a  objects  with 

J  j"1 

repetitions  of  a  choice  permitted.  Thus,  the  number  of  elements  of 
?n+2  corresponding  to  the  composition  (l, l,a^, a^, . . . )  of  n+2  is 


(2) 


1+a2"1 


a0+a  -j  —  1 

2  3 


and  hence 


p(n+2)  = 


al+a2 


a , 


1 


a2+a3~1 


) 


o  «  o 


9 


where  the  sum  extends  over  all  compositions  (a^fa^f  . . . ,a^)  of  n  into 
an  unrestricted  number  of  positive  parts,  and  the  composition  with  one 
part  contributes  1  to  the  sum. 

In  order  to  be  consistent  with  the  defintions  and  notation  of 

Chapter  I  we  write  b(n)  for  p(n+2)  ;  then  the  numbers  b(a,n), 

b  (n),  ...  have  obvious  combinatorial  interpretations.  Using 
Ic 

f (m,n)  =  (m+”‘1)  ,  we  have 
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M 


F(x,y)  =  I  I  (“T1)  An 

t^i  ifci 


xy 

(l"x) (1-x-y)  ' 


and  substituting  this  function  in  (1.15)  gives 


(5) 


B  (x,y)  =  Xy  •  B  (x,  — ^ — )  , 

k+1  1-xy  k  *  1-xy'  9 


for  k  =  1,2,...  .  We  have  B^(x,y)  =  xy/(l-xy)  ;  this  fact  along  with 
(5)  and  an  induction  argument  gives 


(6) 


k 

MX,y)  “  {Ak_1(x)"XyAk_2(x)}  (Ak(x)-xyAk_1(x)  )  ' 


where  AQ(x)  =  A^x)  =  1  ,  and  A^fx)  =  ,  for 

k  =  0,1, ...  .  The  difference  equation  satisfied  by  (A^(x))  can  be 
used  to  prove  that 


(7) 


Xy\  l^x) 

B1(x,y)+...+Bk(x,y)  =  • 


Since  B^+B2+...  =  B  ,  we  put  y  =  1  in  (7)  and  use  the 

difference  equation  satified  by  {A^}  to  find 


(8) 


xA  (x)  \(x) 

B(M)  =  "Vi^  =  ' 1 


Using  the  fact  that  the  limit  of  a  product  is  the  product  of  the  limits, 
(8)  implies 


(4) 
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. 
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{1+B(x, l))^  »  lim 
k  -+oo 


Ak(x) 

\7T^ 


Ak _! (x) 
•^r 


Bf*.1) 

X 


Solving  the  quadratic  equation  in  (9)  for  B(x,l)  and  recalling  that 
b(n)  >  0  ,  we  obtain 


(10) 


*!-•»  '  '  l  *  O  ; 

n 


1.  2n 

thus,  p(n+2)  =  (  ^  )  .  This  generating  function  was  also  found 

in  [6  ]  by  applying  Polya's  fundamental  theorem. 


To  find  b(a,n),  the  number  of  elements  of  pn+2  with  exactly 
a  nodes  in  the  second  level  we  use  (7),  (8),  and  (10)  and  obtain 


(11) 


B(x,y)  =  if111  xy 

k  ->» 


xy ( 1+B(x,  1 )  ] 
l-xy [ 1+B (x, 1) ] 


\  y(l“V l-4"x  ) 

1-  ^  y(lWl~4x  ) 


e 


sl  n 

By  definition  of  B(x,y)  ,  b(a,n)  is  the  coefficient  of  y  x  in 
the  expansion  of  the  last  member  of  (ll)  thus. 


The  number  of  elements  of  P^+2  with  exactly 
the  second  level  and  exactly  k+2  levels  is  bk(a,n)  , 


a  nodes  in 
which  is 
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generated  by  B^(x,y)  given  in  (6).  Similarly,  the  number  of  elements 
^n+2  with  an  arbitrary  number  of  nodes  in  the  second  level  and  exactly 
k+2  levels  is  bk(n)  t  generated  by  Bk(x,l)  . 

The  boundary  conditions  and  difference  equation  satisfied  by 
{A^. }  imply 

hs)  Ak(x)=  (e^+1  -  0g+1)/(e1-e2)  , 

where  (l-0^z) (l-Q^z)  =  1-z+xz^  .  Thus,  for  example,  from  (6)  and  (13) 
we  find  the  explicit  formula 


(14) 


Bfc(X,1)  = 

_ xk(l-4x) _ 

=  (ef +5  +  ef +5  -  xk+1) 


The  first  equality  in  (14)  and  the  difference  equation  satisfied  br 
{A^}  may  be  used  effectively  to  calculate 

2 


(15) 


B1(x,l) 


x 

1-x  9 


B2^X,X^  ~  Tl^)Tl-2x) 


B5(x,1) 


(l-2x)(l-3x+x^) 


B4(x,1) 


4 

x 

( 1 - 3x+x^ ) ( 1 -4x+3x^ ) 


Of  course,  (14)  also  implies  that  each  of  the  sequences  {bk(n)  S  n  =  1,2,...} 
satisfies  a  linear,  homogeneous  difference  equation  of  order  k  with 
constant  coefficients. 
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If  all  of  the  nodes  of  a  tree  have  valence  either  1  or  k+1  , 

then  the  tree  has  kn+2  nodes  altogether.  This  is  because  the  tree  with 

two  nodes  has  this  property  and  every  tree  of  this  kind  with  more  than 

two  nodes  can  be  drawn  by  joining  k  nodes  at  a  time  to  nodes  of  valence 

1c 

1  in  a  tree  already  possessing  this  property.  We  let  ^kn+2  ^enote 
the  subset  of  P^n+g  containing  trees  with  nodes  having  valence  either 
1  or  k+1  ,  and  let  p^(kn+2)  denote  the  number  of  elements  in  this 
set  • 


Each  element  X  of  P,  .  is  such  that  each  level  of  nodes 

kn+2 


except  the  first  level  contains  a  multiple  of  k  nodes.  If  there  are 

ka^  ^  and  ka^  nodes  in  the  (j-l)  and  j  levels  of  X  ,  the 

ttr  s  t 

nodes  in  the  j  level  may  be  joined  to  those  in  the  (j-l)  level 


in  any  of 


(T) 


ways.  This  follows  since  each  block  of  k  nodes 


in  the  j  level  must  be  joined  to  exactly  one  node  in  the  (j-l) 

level,  and  each  node  in  the  (j-l)St  level  is  joined  to  at  most  k 

fch 

nodes  in  the  j  level.  Since  there  are  exactly  k  nodes  in  the 


second  level  of  every  element  of  ^cn+2  *  f°r  n  >  0  ,  the  number  of 
elements  in  ^n+2  *s  b(l>n)  =  p^(kn+2)  y  where 


(16) 


f 


and  the  sum  extends  over  all  compositions  of  n  into  an  unrestricted 


number  of  positive  parts. 


;3...  •  VS d  6:31  ■'  •  -O  '  >.fe  tl 

..  ■  •  ■ 

■■  i  5 
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<3  nJt  adnociald 


•  ('•n-cocos- 


-  67  - 


Using  f(m,n)  =  (k“) 


we  have 


00  00 


(it) 


«”>  •  I  I  O  -V  ■  -  *  , 

Itel  n=l  1**(1+y) 


so  that  (1.17)  becomes 


(18) 


B(x,y)  =  -?&■  +  -L.  r  -*y(l+s )  _  bCx  s) 

'  ’yj  1-xy  +  2jti  /  k  ,,  ,k 

Jc  s  -xy(l+s) 


ds 


It  is  also  easy  to  prove  that 


(19) 


B(x,y) 


ydB(x,0) 

l-ydB(x,0) 


where  d  B(x,0)  denotes  the  i^  partial  derivative  with  respect  to 
y  of  B(x,y)  at  y  =  0  .  This  follows  from  the  fact  that 


(20)  b(i,n)  =  ^  b(l,a1)b(l,a2)...b(l,a1)  , 

where  the  sum  in  (20 )  extends  over  all  compositions  (a^,a^, . . .,a^) 
of  n  into  exactly  i  positive  parts. 

Substituting  the  expression  for  B(x,y)  given  by  (19)  into 
the  right  and  left  members  of  (l8)  and  then  differentiating  with  respect 
to  y  at  y  =  0  in  the  resulting  equation  gives 


(  .  \  i  y  c  *  ■ '  :  H '  ,/\i  VO  •  ■mu- 

l  .  5auf  i  fee  t :  ‘  <  >  :ii  C  »  y  i*  y  oi 
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(21) 


dB(x,0) 


x  + 


1 

2jti 


dB(*.°) 

1-3B(x,0) 


ds 


=  x  {l+6B(x,0)]k  . 


A  purely  combinatorial  proof  of  (21 )  follows  if  it  is  noted 
that  each  element  of  ^kn+2  can  interPreted  as  a  k-tuple 

. . ,X^)  with  X^  e  Pka  +2  9  for  *  =  l>2,...,k  ,  and  where 
(a^,a^, . . • ,a^)  is  a  composition  of  n  into  non-negative  parts. 


Solving  (21)  for  8B(x,0)  involves  the  problem  of  finding 
y  =  8B(x,0)  +1  as  a  power  series  in  x  such  that  xy  -y+1  =  0  \ 
essentially  the  same  problem  was  solved  by  Kelly  [11]  who  used  a 
method  for  inverting  trinomials  due  to  Frame  [2  ].  Following  this 
procedure  we  obtain 


(22)  SB(x,0)=  £  > 

n=l 

so  that  the  number  of  non-map-isomorphic  planted  plane  trees  having 
nodes  with  valence  either  1  or  k+1  is 


(23)  pk(kn+2)  =  (k:1y~i;  (k“)  »  n  =  0,l,  ...  . 

Frame fs  result  is  as  follows;  Suppose  we  have  a  non-associative 
multiplication  defined  on  [a^a^  . . .  ,an)  ;  then  the  number  of  different 
products  which  arise  if  factors  of  a^a  ...an  are  taken  r+1  at  a  time 


1  ii  C  IO l  >5  (IS)  to  io^q  I*lxo**nt<tooi>  yisiwq  A 


,  1  ri.  > 


.' .  .  I  t  I  m.  :q  naa  £>;:)  I.'  r  1  JC-ieso 

■  -  ^  '  :;ri 


■ 
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by  inserting  parentheses  and  keeping  the  order  of  the  a^  unchanged 


is 


kr+1 


(r+l)k 

k 


if  n  =  kr+1  ,  and  is  zero  otherwise. 


From  (lO)  and  (23)  we  see  that  TL  _  =  P0  ^  ,  the  set  of 
non-map-isomorphic  planted  plane  trees  having  nodes  with  valence  either 
1  or  3  i  bas  the  same  number  of  elements  as  ?n+p  •  Sabidussi  [20], 
in  his  review  of  the  paper  of  Harary,  Prins,  and  Tutte  [6  ]  ,  remarks 
that  the  natural  one-one  correspondence  which  they  establish  between  the 
elements  of  these  sets  is  the  most  important  result  in  their  paper.  We 

are  going  to  establish  a  simple  one-one  correspondence  between  the  elements 

2 

of  either  of  these  sets  and  the  elements  of  the  set  _  of  all 

2n-2 

binary  sequences  (a^  a2*  " 9 >a2n  2)  conta*n:*-n8  exactly  n-1  units  such 
that 


(24) 


2j 

\  a.  >  j  ,  for  j  =  1,2,..., n-1  . 

i=l 


The  one-one  correspondence  induced  in  this  way  between  the  elements  of 

T  „  and  P  ^  is  much  simpler  than  the  one  given  in  [6  ]. 

2n+2  n+2 


We  suppose  the  elements  of  PR+p  and  ^2n+2  t0  be  labelled 

in  the  way  described  in  the  third  paragraph  of  this  chapter.  A  unique 

binary  sequence  of  length  2n+3  can  be  assigned  to  each  element  X 

of  P  as  follows:  The  sequence  starts  (1,0,...  )  to  indicate 
n+2 

that  the  root  of  X  has  exactly  one  edge  going  from  it  to  the  node  in 
the  second  row.  Next,  we  add  as  many  l's  to  the  sequence  as  there 
are  edges  going  from  node  i  to  the  nodes  in  the  row  above  it  and 


. 


a  •  -  ■  ^  ■  ->  'roi  o»''i  -  *  -non 


•  *  ‘*9"  '  tvl 

. 


TJtsrh  ni  Hue  j  i  i  in  t:  k  :  s/te  ■  t~*  i  3  o  a  in  awo  I  *? 


j  j  i*  t  b:  -t  .1*  rf3  o  7  rials  io 

:  *)  -.soasups.-  yiJWJkd 


£  >  V  w  ■  >'  *.vriT 
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then  add  0  to  the  sequence,  for  i  =  2,3, . ..,n+2  ,  respectively. 
There  are  as  many  l's  and  o's  in  this  sequence  as  there  are  edges 
and  nodes  in  the  tree  respectively j  hence,  the  sequence  has  (n+l) 
+(n+2)  =  2n+3  digits.  Since,  for  n  >  2  ,  each  sequence  begins 
and  ends  in  the  same  way,  namely  (l,0, 1, . . . ,0,0)  ,  we  can  delete  five 
digits  of  this  sequence  to  obtain  W(x)  =  (a^a^  . .  °,a^2)  .  It  must 
be  that  W(X)  satisfies  (24),  since  (a^a^ . .  is  the  set  of 

digits  which  begins  W(Y)  ,  where  Y  is  a  part  of  X  with  j+2  or 
more  edges.  Furthermore,  if  (a^a^  . .  satisfies  (24)  we  can 

construct  an  element  of  Pn+2  using  (l,0,l,a1, . . ,,a2  ,0,0)  . 


Suppose  X  e  ^2j\+2  anc*  t^iat  X  labelled  in  the  same  way 

as  any  other  elements  of  Pn  _  .  Let  i.  <  i0  <  ...  denote  the  nodes 

2n+2  i  2 

of  X  having  valence  3  and  set  i.  =  j*  •  Corresponding  uniquely  to 

3 

X  is  the  sequence  (<x^,P^)  (a^, P^)  •  •  •  (a-n, Pn)  >  where  the  digits 

and  (3.  correspond  to  the  nodes  joined  to  j*  on  the  left  and  right 
J 

respectively  in  the  level  above  j*  .  A  digit  in  this  sequence  is  0 
or  1  if  the  node  to  which  it  corresponds  has  valence  1  or  3  respeec- 
tively.  It  is  easy  to  see  that  every  element  of  ^2n+2  exacfcly  n 

nodes  of  valence  3  J  also,  it  is  clear  that  n*  is  such  that  a  =  6  =0 
Thus,  we  can  delete  the  last  two  digits  of  (<x^,p^)...  (an>Pn)  without 
loosing  any  information  about  X  ;  hence,  we  define  W*(x)  =  (a^,p^,a2,P2, 

**,an-l,Pn-l^  =  ^al,a2,,,#,a2n-2^  *  Again  w*(x)  must  satisfy  (24) 
because  a  part  of  W*(x)  corresponds  to  a  smaller  part  of  X  having  the 

same  properties.  Furthermore,  it  is  easy  to  check  that  if  (a1,ag, . . «>a2n„2^ 

2 

is  an  elment  of  S„  „  ,  then  we  can  construct  an  element  of  T0  0 

2n-2  dn+d 


53 


. 


t  3  (.i  riaitsa  (,  •  ,,  ,  »)  » 


•  (c.  . .  0  io  •*  ls  r"  ,09 
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using  (ai>a2, . . •^a2n«2,0^°)  •  In  Figure  1,  the  correspondence  between 

2 

the  elements  of  S  0,  P  0  ,  and  T  is  shown  for  n  =  3  . 

2n-27  n+2  2n+2 


(1,0,0, 1) 


r 


V4 


2 

10 


(0,1, 1,0) 

5 


'V 


3 
2  0 
lO 


0 


Q 


2 

16 


Figure  1, 


The  correspondence  between  the  elements  of  S 


2n-2 


and  T 


2  n+2 


is  a  special  case  of  a  more  general  result.  Let  ^  denote  the  set 

of  binary  sequences  (a^a^,  . .  .,8^  containing  exactly  n-1  units 

such  that 


(25) 


kj 

z 

i=l 


ai  >  J  , 


for 


j  =  1,2, . . .,n-l 


There  exists  a  natural  one-one  correspondence  between  the  elements  of 
k  k 

^kn  k  an<*  ^kn+2  °  '^ie  Pro°f  °f  this  fact  is  a  straightforward 

generalization  of  the  argument  presented  in  the  last  paragraph.  This 
correspondence  for  k  =  3  and  n  =  2  is  illustrated  in  Figure  2. 
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Figure  2. 


We  give  the  name  cypress  tree  to  those  elements  of  Pr  having 
the  property  that  the  first  node  on  the  left  in  each  level  is  joined  to 
the  first  node  on  the  left  in  the  level  above  it:  also,  C  denotes 
the  set  of  all  cypress  trees  with  n  nodes.  By  the  same  kind  of 
argument  used  to  prove  (3),  it  follows  that  c(n+2)  ,  the  number  of 
cypress  trees  with  n+2  nodes, is  equal  to 


(26) 


b(a) -  y  ( vv2 ^ )...( 3i-i+3i 

1— i  \  ap"l  '  a^-l  /  \  a. -1 


W* 


a .  -  +a . -2 


where  the  sum  extends  over  all  compositions  (a^,a^, . . «,a^)  of  n  into 
an  unrestricted  number  of  positive  parts. 

Since  f(m,n)  =  >  we  have 


(27) 


F(x,y) 


00  00 


m=l  n=J 


/m+n-2  \ 
1  n-1  } 


m  n 
x  y 


22L 


1-x-y 


hence,  substituting  this  into  (1,15)  gives 


.£  3  .1.  ;  :  '•  ' 


lo  i  li  aor sXa  »ao rfj  03  >93:3  .  -v  *r  yo  smarj  >rf-j  sv  :g  sW 

f  d:  o  i  » I  1/5  ,  •  t>iij  3  c  )  y:  qoxq  »rii 
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N 

eon  +  ,;i73  avaiqyo 


,:? )  l  •  <■ 


(as) 


avari  aw  ,  (  ,  )  «  (nta)^  aoal8 

-  <r,*)t  (TS) 
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(28) 


Bk+l(x,y)  =  xyBk(x,  . 


We  have  B^(x,y)  =  xy/(l-xy)  ,  and  a  simple  induction  using  (28)  gives 


(29) 


B.  (x,y)  = 


k 

X.Z 


Ak(x)-xyAk^1(x) 


where  the  A^'s  are  the  polynomials  defined  after  (6).  From  (29)  it 
is  evident  that  {t>k(n)}  ,  where  b^(n)  denotes  the  number  of  cypress 
trees  with  exactly  k+2  levels  and  n+2  nodes,  is  generated  by 
Bk(x,  1)  =  x^A^+^(x)  ,  a  rational  function.  For  example,  when 
k  »  1,2,3>^>  and  5  >  B^(x,  l)  is  respectively 


(50) 


x 


X 


1-x  9  l-2x  9  2  9  ,  2  9 

l-3x+x  l-4x+3x 


and 


l-5x+6x^-x^ 


The  third  of  these  relations  implies  b^(n)  =  F^n  where 
(F^)  denotes  the  Fibonacci  sequence;  thus,  the  number  of  cypress  trees 
with  exactly  five  levels  and  n+2  nodes  is  F~,  for  n  =  1,2,...  . 

2 

Suppose  a  and  3  are  real  numbers  such  that  t  -at +3  =  0 

has  complex  roots;  then  it  is  known  that  the  Taylor  expansion  of 
2  -1 

(l-az+3z  )  has  infinitely  many  negative  coefficients.  Now  since 

00 

(51)  — — 2  =  y  > 

1-Z+XZ  l—i 

k=o 

we  must  have  Ak(x)  <  0  for  infinitely  many  k  whenever  x  =  l/(4-€)  , 
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4  >  €  >  0  .  But  x[?]  A^(^)  =  +  ...  ,  so  for  some  k  , 

)  =  0  has  a  real  root  arbitrarily  close  to,  but  less  than  4  . 
This  imples  b^Cn)  >  (^-€)n  ,  for  sufficiently  large  k  and  n  . 


tl  <  jf»  .  C  '  '  •  ,/ 
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CHAPTER  VI 

SOME  APPLICATIONS  TO  THE  THEORY  OF  NUMBERS 

It  is  sometimes  of  interest  in  the  Theory  of  Numbers  to 
evaluate  or  estimate  sums  which  extend  over  all  or  only  a  special 
subset  of  the  compositions  of  n  .  Perhaps  the  most  common  sum  of 
this  type  is 


(1)  r (n)  =  V  s (d)  , 

ajn 

where  the  sum  extends  over  all  compositions  of  n  with  the  parts 
equals  that  is,  the  sum  extends  over  all  positive  divisors  of  n  . 
The  formulas  of  Chapter  I  can  be  specialized  in  order  to  deduce  some 
of  the  well  known  facts  about  generating  series  of  sequences  (r(n)) 
defined  as  in  (l). 


Suppose  W(x) 
respectively;  then  define 

(2)  f(m,n)  =  • 


and  F(x) 
F(x,y) 

r  f(n), 


0 


generate  (w(n))  and  {f(n)) 
and  G(x)  in  (1.7)  by 


if  m  =  n 
otherwise 


,  and 


(3)  g(n)  =  f (n)  w(n)  . 

For  any  composition  . . ,,a^)  of  n  we  have 


&i  •&  V'  wm:*  i  1  ?  '• ,  ;'  5  '  '  :i 

:  V  '  '  "  •  -■  ■'  '•  I'- 
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1  .  •.  -  «>  .  •  I  s  -  *  '■  •  '■  ,l 


^  X 

i  £  ‘XM  I  3  -J8  -  ^  Siriw 
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(a)fn/a(a),  if  a |n 


w 

f(a.,a  )  f(a  ,a  )  ...  f(a  ,a  )  g(a  )  =J 

\  /  7  | 

so  that 

J  L  0  , 

otherwise 

w(a)fn/a(a),  if  a  In 


(5) 

b(a,n)  = 

0  ,  otherwise 

,  and 

(6)  b (n)  =  V  w(d)fn/d(d)  . 

djn 

The  definitions  in  (2)  and  (3)  also  imply  F(x,y)  =  F(xy)  , 

and 

00 

(7)  G(x)  =  f(n)w(n)xn  =  7^-  J  F(xs)w(~)  , 

n=l 

where  c  is  a  contour  in  the  s  plane  which  includes  the  singularities 
of  W(^-)/s  ,  but  excludes  those  of  F(xs)  . 

The  case  when  f(n)  =  1  ,  for  n  =  1,2,...,  is  of  particular 
interest.  Here  F(x)  =  x/(l=x)  and  G(x)  =  W(x)  ;  thus,  (1.13)  implies 
B^(x,y)  =  W(xy)  ,  and  (1.15)  implies 


Bk+i(x>y) 


xyBk(x,s) 

—^7 - -t—  ds  =  B,  (x,xy) 

s(s-xy)  kv  ' 


From  (8)  we  deduce  that  Bk(x,y)  =  W(x  y)  ,  and  hence  that 


B(x,y)  =  W(xy)  +  W(x2y)  +  W(x^y)  + 


(9) 
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Under  present  assumptions,  the  functions 


BJ(x,y) 


are  easily 


found  since  the  system  (1,23)  becomes 


(10) 


i  p  dPBj  p 

—  =  w(p)*p  +  —  X 


for  p  =  l,2,...,j  ,  so  that 


(u) 


dpB^  _  w(p)xp 

P!  =  (X-xP) 


From  (11)  we  deduce  that 


(12) 


BJ(x,y)=  V  , 

P4i  u-*p) 


and  letting  j  tend  to  infinity  in  (12)  we  obtain 


(13) 

00 

B(x,y)  =  V 

»(p)*pyp 

p4i 

(l-xP) 

written 

Combining  (5),  (9), 

00 

and  (15),  the  final  result  may  be 

(14) 

B(x,y)  =  y 

V  w(d)ydxn 

n=l 

djn 

00 

00 

.  y 

w(p)xPyP  =  V  W(xny)  . 

pi 

(1*xP)  n4 

Umco  »»  (St)  l  x  •  fc,‘  '  51,1  00“ 


n9l“’W 

WO 
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From  (14)  we  can  deduce  many  well  known  results  (see  pages 

257“258  of  Hardy  and  Wright  [9])»  For  example,  if  w(n)  =  p(n)  , 

the  Euler  ^-function,  it  is  known  that  \  p(d)  =  n  ;  hence,  (14)  implies 

din 


(15) 


00 


00 


4  i  -  1 

n=l  d|n  n=l  v  7 


,  and 


(16) 


00 


00 


(l-x)‘ 


n=l  dYn  n=l  u  } 


that 


Also,  if  w(n)  =  n(n)  ,  the  Mobius  function,  it  is  known 


(IT) 


y  i*(d)  -  { 1  ’ 

dYn  U’ 


l  ,  if  n  >  1 
otherwise 


thus,  (14)  implies 


00 


00 


(18) 


I 


,  \  n  n 

p(n/.x.  y 


p,(d)y^xn  =  \  J  and 

-1  HI  4  (1-Xn) 

n=l  d  In  n=l  v 


00 


00 


(19) 


X  = 


n=l  d|n  n=l 


i  r  i 

Finally,  if  w(n)  =  n  ,  \  d  =  crk(n)  ;  usually,  we  write 

dTn 


-•/  /  ,  '  •'  1 


' 


(*')  \  ■'  1  :  5  V 


'-N 


K1& 


' 
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c0  =  t  ,  cr  =  o-  .  For 
is  x(l-x)  *  ,  x(l-x)  ^ 
(14)  implies 


k  =  0,1,2,  the  generating  function  of  {n  } 
,  and  (x+x  )(l-x)  respectively;  hence, 


00  00 


00  00  00 


(22) 


00 


00 


Zr — I  2  n  r— I  n  £ 

a9(n)x“  =  V  -2-2-  =  V  2^2_ 

^  A  (l-xn)  Zj  fl-xn 

n=l  n=l  U  X  1  rfel  X 


00 


n  2n 
:  - 

(!-*“)• 


Another  class  of  problems  in  the  theory  of  numbers  asks  for 
the  number  of  compositions  of  n  which  satisfy  certain  requirements. 
Some  of  these  problems  stipulate  how  many  parts  may  be  used  and  which 
numbers  may  serve  as  parts  of  the  composition;  a  generating  function 
for  compositions  of  this  kind  is  easily  obtained.  Let  b(n)  denote 
the  number  of  compositions  of  n  which  satisfy  these  stipulations 
and  define 


(23) 


p(n) 


if  a  composition  may  have  n  parts 

) 

otherwise 


if  n  may  serve  as  a  part  in  a  composition 


t 


w(n)  = 


0 


otherwise 


..  .  £ •<  1  (  I) 


(Vi 


f  »T 


00 

t(f.  )T> 

f:n 


;;  j  '  «*(,  )c»  ;  is*> 

C(Vx)  c  -  J>n  |An 

*  *  >f  .  r.-„  r  3  1  -s  H  m  f.  ?q  ^  8  -  ’ :A 

,  ,  ■  ■.  L.,v.  .  .  >■ ■<■■■-.,  wo,  ^  *  5  08 

■ 

f  ..-a  \-u .  ■  ‘  ••  «1  t:  fH'  5  ■  7  ■  !i  q'  '  IOi 

$nJtl»b  bn* 
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Suppose  P(x)  ,  W(x)  ,  and  B(x)  generate  {p(n)]  ,  {w(n))  ,  and 

(b(n) )  respectively,  then  it  is  easy  to  show  that  B(x)  =  P(w(x)) 
(for  a  proof,  see  Klarner  [13];  this  result  is  a  generalization  of  a 
theorem  proved  by  Moser  and  Whitney  [14]). 

Problems  involving  enumeration  of  compositions  become  more 
difficult  when  in  addition  to  stipulating  the  number  and  type  of 
parts,  the  pairs  of  consecutive  parts  are  restricted  in  some  way  as 
well.  A  common  type  of  composition  satisfying  properties  of  this 
kind  is  a  partition  of  n  which  we  will  consider  presently. 


Suppose  a  set  S  of  positive  parts  and  a  subset  P  of 
S  X  S  are  given  and  define 


(24) 


f (m,n) 


if  (m, n)  e  P 
otherwise 


(25) 


s(n) 


r  w(n)  ,  if  n  e  S 
l  0  ,  otherwise 


where  w(l),  w(2),  ...  are  the  weights  of  1,2, ... 


Now  defining  b(n)  in  the  usual  way  we  have 


(26) 


b(n)  = 


f  (ax 


>  ) 


f(a2>a5) 


£^i.vai)  8(ai} 


^  w(ax)  w(a2)  . . .  w(a1) 


f 


t  ^)w)  t  99  (*)S  c-  a  ,  x  H  t  ^  ^  s»aoqqu3 


•-«;  »  (x)  *e;tt  *v  Of:  ’  .c-j  ai  32  v23.o  (  [a)  1) 


‘  t  * ,  '  •  .  . 


9.  rooad  ciiv  .  1  ;oo  e  o  :t  .  /.i» 

os-  -  ir  m  gnidii  -  •  >rtv  Oluoil  i!_. 

;  -  s  •: 

e.-.rfJ  :o  m  l:  -qr  .  i.  .1  11  .  .  w  00  *  o  no  \>>  1  , XI  n 
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where  the  first  sum  extends  over  all  compositions  of  n  and  the  second 
sum  extends  only  over  those  compositions  (a^,a^, . . .,a^)  of  n  such 
that  a.  e  S  and  (a  ^,a.)  e  P  ,  for  j  =  l,2,...,i  .  Thus,  the 
general  problem  we  have  set  out  to  solve  may  be  treated  using  the  methods 
of  Chapter  I. 


For  example;  to  find  the  generating  function  of  {b(n)}  , 
where  b(n)  denotes  the  number  of  ways  of  expressing  n  as  a  sum  of 
one  or  more  consecutive  numbers,  we  put  g(n)  =  1  ,  and 


(27) 

so  that 


f(m,n) 

G(x)  =  x/(l-x)  , 


F(x,y) 


if  m  =  n+1 

t 

otherwise 

2 

x  y/(l-xy)  and 


(28) 


b(a,n) 


1  ,  if  n-a  =  l+2+...+i  ,  for  i  >  0 
0  ,  otherwise 


Using  (1.15)  we  get 


(29) 


Bk+1(x’y)  -  m  L 

J  c 


2  2 

x  y  B,  (x,s) 


k 

s (s-xy ) 


ds  =  xyBk(x,xy)  , 


so 


that  evidently  B,  (x,y)  =  y^/(l-x^y)  ,  and  hence 


Y  xk(k+l)/2  yk 

k4 


(30) 


B(x,y)  = 


*  .  .  7  r.s  2<  :v;  .  v  )  . 

»»•  I  ,  us  l  -  iff  .  I  t  ..a  3f,rfJ 

* 

t  I  (■•?  ';<  /  q  .  7  ...vi  .i'u.tr;  O'  •  I  •  7  0  9H0 


J  x-  -  ii 


9a.iwTl-7rjO 


t  0-J 


X 

/'•  ,£  (vx  I !  \Y  (  .  ■  V.  -  -ri-7  V3 


*  +  .  ,  -  :  1  8  i  f!  . 

V'  *.LjrJ3o 


E)d 


-  >  \  l  ' 


(P£: 


rabivo  lari* 
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To  find  the  generating  function  of  (b(n))  ,  where  b(n) 
denotes  the  number  of  compositions  of  n  with  contiguous  parts  distinct, 
we  set  g(n)  =  1  ,  and  define 


(31) 


f  (m,  n) 


if  m  ^  n 
otherwise 


in  this  case  we  have  G(x)  =  x/(l-x)  and  F(x,y)  =  xy (x+y-2xy)/(l-x) (l-y) (l-xy)  . 
Substituting  these  functions  into  (1.17),  integrating,  and  transposing 
a  term  gives 


(32) 


B(x,xy)  +  B(x,y)  =  {1+B(x,l)j  ^ 


EL 


xy 


£ 

Now  we  write  x  y  for  y  in  (32),  multiply  the  resulting 
identity  by  (-1)  ,  for  k  =  0,1,..., K  ,  and  sum  over  k  to  obtain 


k  k+1 


(33)  B(x,y)  +  (-l)KB(x,xK+1y)  =  (1+B(x,l)}  V  ^ 

to  Ux  y 


K+1 

But  B(x,x  y)  tends  to  zero  as  K  tends  to  infinity  for 
x|  sufficiently  small;  thus  (33)  implies 


(34) 


B(*.l). 

1+B(x,l) 


00 


(  i \^+l  k 

k  *  =  T(x)  • 

1-x 


Using  (23),  we  see  that  B/(l+B)  generates 


(35) 


V  (-1)1  b(a1)  b(a2)  ...  b^)  , 


(n)d  9 tai:  [(a)d’.  Jo  aoUiOiA  1.  »»»•;  «  iwii  .o 

' 

,  ‘ 
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M  > 


' **- '  J-. . 
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where  the  sum  extends  over  all  compositions  of  n  \  also,  (lb) 
implies  T(x)  generates 

(36)  V  (-l)d+1  =  TQ(n)  -  Te(n)  =  T*(n)  , 

djn 

where  T0(n)  and  Te(n)  denote  the  number  of  odd  divisors  and  the 
number  of  even  divisors  of  n  respectively.  Thus,  if  2k  is  the 
highest  power  of  2  which  divides  n  ,  T*(n)  *  (l-k)  T(n/2k)  . 

The  relation  (34)  also  implies  B(x,l)  =  T(x)/(l-T(x)  ,  so 
that  using  (23)  again,  we  have 


(37) 


T*(ax)  T*(a2)  ...  T*(at) 


i 


where  the  sum  extends  over  all  compositions  of  n  .  Equating  (35)  and 
(36)  we  have  the  inverted  form  of  (37) s 


(38) 


T*(n)  = 


^(-l)1  b(fll)  b(a2)  ...  b(ai)  , 


where  the  sum  extends  over  all  compositions  of  n  .  The  sequence  (b(n)} 
has  in  part  the  form  [1,1,3,4,7,14,...)  . 

The  number  of  compositions  of  n  in  which  no  more  than  one 
part  of  each  pair  of  consecutive  parts  is  1  can  be  found  by  setting 
g(n)  =  1  and 


f  (m,  n 


0  , 
1 


if  m  =  n 


1 


(39) 


f 


otherwise 


r)  ,  idii 


;>  • 


(<*) 


jIoijo'-j  (n)  i 

•v 
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Using  G(x)  =  x/(l-x)  and  F(x,y)  =  (x+y-xy)/(l-x) (l-y)  in  (I.I7), 
we  have 

(^°)  B(x,y)  =  jSL-  {1+b(x, l) )  -  xy  dB(x,0)  , 

where  dB(x,0)  denotes  dB(x,s)/ds  at  s  =  0  .  Differentiating  with 
respect  to  y  at  y  =  0  in  (40),  or  putting  y  =  1  in  (40)  gives  two 
equations  linear  in  B(x, l)  and  dB(x,0)  ;  solving  for  B(x,l)  gives 

(41)  B(x, l)  =  — ^  =  x+x2+3x\5x^+9x^+  ... 

1-x-x  -yp 

It  is  interesting  to  note  the  similarity  between  this 
generating  function  and  the  generating  function  given  by  taking  s  =  1 
in  (4.13)  which  generates  the  number  of  two  row  n-ominoes. 

Now  we  are  going  to  show  how  the  common  identities  involving 
the  generating  functions  of  various  partitions  of  n  can  be  made  to 
follow  from  the  formulas  of  Chapter  I.  It  has  already  been  observed 
in  Chapter  IV  that  another  parameter  can  be  introduced  which  keeps 
track  of  the  contribution  to  the  sum  which  defines  b(n)  made  by 
compositions  of  n  with  exactly  i  parts.  To  do  this,  we  simply 
write  zF(x,y)  and  zG(x)  in  place  of  F(x,y)  and  G(x)  ;  then, 
for  example,  (1.15)  becomes 

(42)  B  (x,y,z)  =  2 ^  f  F(xy,i)  Bk(x,s,z)  *2. 

J  c 
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■  ■  ■ 
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and  hence. 


(^3)  B(x,y,z)  =  B1(x,y)z  +  B2(x,y)z2  +  ...  , 

where  B^(x,y)  is  defined  in  (1.15).  Also,  it  follows  from  that 


(MO 


B(x,y,z)  =  zG(xy)  +  f  F(xy,j)  B(x,y,z)  —■ 

•'  c 


We  note  for  purposes  of  comparison  that  the  following 

technique  gives  a  generating  function  in  a  certain  form  for  most  kinds 

of  partition  functions.  Suppose  0  <  a^  <  a2  <  . . .  is  a  given  set 

of  integers  and  let  [a±i’Cii2’  *  *  * )  be  a  set  non-negative  integers 

associated  with  a^  ,  for  i  =  1,2,...  .  The  number  of  partitions  of 

n  involving  exactly  k  of  the  parts  a^,a^,...,  where  a^  is 

selected  as  a  part  exactly  times,  for  some  j  and  i  =  1,2,...  , 

k  n 

is  the  coefficient  of  z  x  in  the  expansion  of 


(M>) 


x^  } 


a, 


'2v 


a  a„ 

2  2v 


}••• 


For  examp le ,  if  {cl ^ ^ , ol^2, *  *  *  3  ~  {0,1,...}  y  for  i  —  1,2,..., 
the  generating  function  is 


(MS) 

IT  fl+ZX  i+ZdX  ■*■+...)  = 

IT  1 

‘  i  .  a .  ,  * 

i=l  i 

L=1  (l-zx  L) 

also,  if 

[a^^,a^2, . . . }  =  {0,1}  ,  for 

i  =  1,2,...  ,  the  generating 

function 

is 

the  generating 
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(V7) 


It  follows  from  (46)  that  pk(n)  ,  the  number  of  partitions 
of  n  into  exactly  k  positive  parts,  is  generated  by 


(48) 


P(x,z) 


1 

(l'-zx^) 


If  g(n)  =  1  ,  and 


m 


f  (m,n) 


r  1  ,■  if  m  <  n 
4  0  ,  otherwise 


f 


then  b(a,n)  is  the  number  of  partitions  of  n  with  smallest  part 
equal  to  a  .  Alternatively,  if  g(n)  =  1  ,  and 


(50) 


£(m,n) 


r  1  ,  if  m  >  n 
4  0  ,  otherwise 


then  b(a,n)  is  the  number  of  partitions  of  n  with  largest  part 
equal  to  a  . 


Using  the  definition  of  (f(m,n))  given  in  (49),  we  have 
F(x,y)  =  xy/(l-y) (l-xy)  ;  substituting  this  into  (1.15)  we  obtain 
a  recurrence  satisfied  by  (Bk(x,y))  which  along  with  BL(x,y)  = 
xy/(l-xy)  implies 


Bk(x,y) 


k 

x  y  _ _  . 

(l-x)  ...  (1-Xk"’^)(l-Xky) 


(51) 


t  'v  » V  '  ^  ^  ^ 
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hence, 


00 


(52) 


B(x,y,z)  =  V  — 
A  (l 


k  k 
yx  z 


^  (l°x)  o.o  (l-xk  1)(l-xky) 


Similarly,  the  definition  of  {f(m,n))  given  in  (50)  implies 
F(x,y)  =  xy/(l-x)(l-xy)  ,  and  this  leads  to 


(53) 


hence, 


k 


B  (x,y)  =  - —  o  .  , 

(l-xy)(l-x  y)  .0.  (l-xky) 


00 


(5*0 


B(x,y,z) 


=  y  - _ ^ 

M  ( 1 -xy ) ( 1 


k  k 
z 


^rj  (i-xy)(i-x2y)  (i-xky) 


When  y  =  1  in  (52)  or  (54)  the  generating  functions  are 
both  equal  to  P(x,z)  -  1  ;  hence  we  obtain  Euler's  famous  identity 


00 


(55) 


00 


n 


i-l  (l-zx1) 


=  1  + 


k  k 
x  z 


(l-x)(l-x2)  ...  (1-Xk) 


If  we  want  the  generating  function  for  the  number  of  partitions 
of  n  with  the  smallest  part  restricted  in  some  way,  we  modify  (51) 
and  (52)  by  redefining  G(x)  .  The  generating  functions  (53)  and  (54) 
can  be  modified  similarly  to  treat  restrictions  on  the  largest  part  in 
the  partition. 

Of  course,  [f(m, n)}  and  [g(n)j  may  be  defined  in  many 
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ways  so  that  b(n)  enumerates  a  variety  of  partitions  of  n  .  Using 
the  technique  of  Chapter  I  ,  it  generally  develops  that  the  functions 
B^(x,y,z)  tend  to  a  generating  function  having  a  form  given  by  (45), 
while  the  function  B(x,y,z)  =  B^(x,y)z  +  B^(x,y)z  +  ...  will  usually 
have  a  form  similar  to  (52)  or  (54).  Equating  these  results  sometimes 
leads  to  a  classical  identity  due  to  Euler,  Cauchy,  Jacobi,  et  al.  ; 
generally,  the  methods  used  in  the  first  place  in  establishing  these 
identities  are  simpler  and  more  elementary  than  an  application  of  the 
methods  of  Chapter  I.  We  have  only  mentioned  this  application  to 
indicate  a  connection  between  the  generating  functions  of  partition 
functions  and  the  Fredholm  equation  which  is  the  connecting  link 
between  all  of  the  problems  discussed  so  far. 
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CHAPTER  VII 

GENERALIZED  GRAPHS 

Let  P^(x)  denote  the  collection  of  ordered  k-subsets  and 
pk(X)  the  collection  of  ordered  k-selections  of  the  set  X  .  Let  G 
be  a  permutation  group  of  degree  k  whose  elements  permute  the  order  of 
the  terms  in  any  k-tuple.  Two  k-tuples  are  G-equivalent  if  there  is  a 
permutation  in  G  which  transforms  one  into  the  others  clearly,  G 
partitions  P^(X)  an£*  ^  (X)  into  collections  of  equivalence  classes 
which  we  denote  by  GP^(X)  and  GP^(X)  respectively.  For  example,  if 
G  =  ,  the  symmetric  group  on  k  letters,  S^P^(X)  denotes  the 

collection  of  all  unordered  k-subsets  of  X  since  two  k-tuples  are  the 
same  if  and  only  if  they  contain  the  same  elements. 

Now  let  G  be  a  given  subgroup  of  S,  and  let  N  =  [l,2,...,m]  ) 

lc  m 

we  call  the  elements  of  N  nodes  and  the  elements  of  GP,  (N  )  edges. 

m  .  — k  m  v 

We  write  G  (m)  for  the  set  of  all  (k,G)~graphs  on  m  nodes  where  a 

it 

(k,G)-graph  consists  of  the  set  of  nodes  N^  and  some  subset  E  of 
GPk^Nm^  called  the  edge  set  of  the  graph.  The  subset  of  Gk(m) 
containing  graphs  with  exactly  n  edges  will  be  denoted  by  G^(m,n)  . 

* 

If  k  is  a  permutation  of  N  ,  it  induces  a  permutation  it 

of  GP,  (N  )  defined  by 
k^  nr 

Jt  4  (a|f  » » 3^)  “*(jt(a^)j°"°>Jt(ak))  f  3  ^k^m^  4 
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Two  graphs  in  G^(m)  are  isomorphic  if  there  is  a  permutation  of  the 
nodes  of  one  which  transforms  its  edge  set  into  that  of  the  second 
graph.  The  set  of  non- isomorphic  (k,G)-graphs  with  m  nodes  and  n 
edges  will  be  denoted  by  Gf  (m,n)  and  we  define  G*(m)  =  U  G*(m,n)  . 

K  tC  n  K 

Let  E  denote  the  edge  set  of  a  graph  A  €  G^(m, n)  and 

define  A  e  Gm_^(m, n)  as  the  graph  having  the  edge  set  E  =  {Nm»x|x  e  E] 

This  one  to  one  correspondence  shows  that  |Gk(m,n)|  =  |G^  k(m,n) |  ; 

hence,  |G^(m)  |  =  |G^  jc(m)  |  •  Next,  if  A,  B  €  G^(m,n)  are  isomorphic, 

then  A  and  B  are  isomorphic  as  well.  To  see  this,  suppose  jt  is 

a  permutation  of  N  such  that  j?  transforms  the  edge  set  of  A 

m  0 

# 

into  that  of  B  ;  then  at  x  =  y,  x  e  A  ,  yeB  implies 
* 

rt  =  Proves  |G*(m,n)|  =  jG*  ^(m,n)j  and  from  this 

it  follows  that  |G*(m) j  =  |G*  ^(m) |  • 

When  k  =  2  ,  G  may  be  chosen  as  either  of  two  groups, 
the  identity  group  I  of  degree  2  or  the  symmetric  group  S  of  degree 
2  .  The  edges  in  IP2(m)  and  SP^m)  are  said  to  be  directed  and 
undirected  edges  respectively;  also,  (m )  ,  S^m)  ,  I*(m)  ,  S*(m) 

are  called  the  labelled  directed,  and  linear  graphs  on  m  nodes,  and 
the  non-isomorphic  directed,  and  linear  graphs  on  m  nodes,  respectively. 
We  are  going  to  consider  the  problem  of  enumerating  the  sets  I^(m)  , 

S,  (m)  ,  I*(m)  ,  and  S*(m)  ,  where  I  and  S  denote  the  identity 

K  1C 

and  symmetric  groups  of  degree  k  ,  respectively.  To  do  this  we  will 
require  Polya's  enumeration  theorem  which  we  state  here  for  convenience. 
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Suppose  S  is  a  set  of  distinct  objects  and  that  w(n) 
objects  have  been  assigned  the  weight  n  ,  for  n  =  0,1,...;  we 


define  f(x)  =  \  w(n) 


n 


x 


If  G  is  a  permutation  group  of  degree  k  , 


£ 

then  the  weight  of  an  equivalence  class  in  GP  (s)  is  the  sum  of  the 
weights  of  the  objects  contained  in  a  representitive  k-tuple  in  the 
equivalence  class.  Letting  W(n)  denote  the  number  of  equivalence 
classes  with  weight  n  ,  for  n  =  0,1,...  ,  we  define  F(x)  =  ^Tw(n)xn 

P61ya’s  theorem  states  that 


(2) 

F(x)  =  Z(G:f(x),  f(x2), 

f(xk))  , 

where 

(3) 

z(Gsf1,f2,0°0,fk)  =  "|g[ 

y  =(#)  , 

jteG 

and  c(jt) 

•  •  « 

Jj  J o  Jk 

=  f .  f.  f.  if  ye  has 

12  k 

j.  cycles  of  length 

i  =  1,2,. 

ao,k  ;  Z(G; f . . o , fk)  is  the 

cycle  index  of  G  . 

i  for 


Before  proceeding,  we  note  that  if  the  number  of  edges  in 

GP,  (N  )  is  known,  then  it  is  easy  to  find  the  number  of  graphs  in 
k  m 

Gk(ra)  and  Gk(m,n)  .  In  fact,  if  lGpk(\)  I  *  N  >  we  have 
|Gk(m,n)|  =  (^)  and  |Gk(m)  |  =  2N  .  Thus,  since  |IPk(Nm)|  =  kUk) 
and  |SPk(Nm)|=  (“)  we  have 
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For  the  time  being  we  assume  that  G  is  the  identity 


group  I  of  degree  k  )  thus,  IP^N  )  denotes  the  set  of  ordered 
subsets  of  N  .  Wfe  want  to  use  P6lya's  theorem  to  determine  the 
number  of  elements  in  I*(m, n)  and  I*(m)  .  In  the  application  of 
the  theorem,  the  set  S  consists  of  the  two  objects  "is  not  contained 
in  the  edge  set"  and  ”is  contained  in  the  edge  set”  which  we  weight 
0  and  1  respectively.  If  the  elements  of  IP^N^)  are  put  in 


some  fixed  order,  say 


IP 


k!(“) 


,  ,  ,,,,X 

'  2  k!(m 

n 


each  element  of 


(s)  may  be  interpreted  as  a  description  of  a  graph  contained 

th  k!(£) 

in  by  saying  the  i  component  of  X  e  IP  (s)  tells 


whether  X,  is  or  is  not  contained  in  the  edge  set  of  the  (k, l)- 
graph  being  described.  The  permutation  group  G  of  IP  (S)  is 
induced  by  the  symmetric  group  on  as  described  in  (l).  To  apply 

the  theorem,  we  need  the  cycle  index  of  G  and  use  f(x)  =  1  +  x  . 


Let  denote  a  permutation  of  with  cycles  of 

length  i  involving  the  ij^  nodes  contained  in  for  i  =  1,2, ...,m  , 

* 

We  want  to  determine  the  number  of  cycles  of  length  i  in  it  for 

/in  \  ^ 

i  =  1,2, ...,k!(  )  .  Note  that  it  maps  an  edge  having  a.  nodes 
in  onto  another  edge  having  a^  nodes  in  for  i  =  l,2,...,m  . 

It  is  easy  to  see  that  the  edge  cycles  induced  by  it  among  the  set 
of  edges  having  a  nodes  in  J„  with  a  +  a  +  ...  +  a  =  k  , 

1  1  I  d. 
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>  0  ,  have  length  A(a^)  ,  the  least  common  multiple  of  the  sub¬ 
scripts  of  the  non-zero  a's  .  Furthermore,  if  the  non-zero  a's 

are  a.  ,  a.  , . „ . ,  then  the  number  of  edges  of  this  kind  is 
1  l2 


(6) 


k! 


Vi. 


a 


Vi. 


=  v^)  . 


This  follows  from  the  fact  that  there  are 


iJ. 


a 


ways  to  select 


a^  elements  from  for  i  =  and  k!  ways  to  distribute 

these  elements  in  a  k-tuple.  Thus,  «  induces  v(a^)/A(a_^)  edge 
cycles  of  length  A(a^)  in  the  set  of  edges  having  a^  nodes  in 
for  i  =  l,2,.o.,m  in  each  edge0  Letting  (a^, a^, . . . , a^)  range  over 

all  compositions  of  k  into  non-negative  parts,  we  obtain  a  complete 

* 

list  of  the  cycles  of  it 


jf  ^2  *  ^2 

If  we  write  c(  jt  )  =  f^  f  =  o.  and  c(  jt  )  =  gg  ••• 

explicit  formulas  can  be  obtained  using  the  result  of  the  last  paragraph. 

*  , 

For  example,  the  contributions  to  c(  )  when  =  k  ,  or 

a  +a  =  k,  a,  a  /  0  ,  or  a  +  a  +a^  =  k,  a  ,  a  ,  a.  ^  0  , 
r  s  r  s  '  r  s  t  7  r  s  t  ' 

are  respectively 


j  .•  :  -  3  '  •  . 


- 


•  .  '  e’< 


' 


-  9^  - 


(8) 


f  r  f  S 


‘L( 


SJ. 


a 


s] 


g 


[r,s] 


(9) 


jr  js  jt 
f  r  f  s  f  c 
r  s  t 


“I  C:X  roc1')/— 


g 


[r,s,t] 


r,  s 


where  the  sums  in  (8)  and  (9)  extend  over  all  compositions  (a^  a  ) 

and  (a  , a  ,  a  )  of  k  into  positive  parts, 
r  s  t 


Using  these  formulas  it  is  not  difficult  to  find  the  cycle 
index  involved  when  m  =  6  and  k  =  3  ;  it  is 


(10) 


{g}2°  +  60g2°  +  40g^°+  l8og^°  +  Ikbg^  +  I20g 

+  +  1+,°8183^  +  120g283865  J  * 


20 

6 


The  number  of  non- isomorphic  (3*l)“graphs  with  n  edges  and  6  nodes 

n  £ 

is  the  coefficient  of  x  in  (10)  when  1  +  x  is  substituted  for 


8i 


Now  we  assume  G  =  S  ,  the  symmetric  group  of  degree  k  ; 
thus,  elements  of  SP,  (N  )  are  simply  k-subsets  of  N  .  It  is 


kv  m 

easy  to  see  that  |S*(m,n)|  =  1  and  |S*(m) |  =  m+1  .  The  generating 
series  for  |S*(m,n)|  ,  the  number  of  linear  graphs  with  n  edges, 
was  given  by  P6lya  [16]  ]  subsequent  expositions  of  this  result  have 
been  given  by  Harary  [8]  and  Riordan  [19].  The  cycle  index  of  the 
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permutation  group  induced  in  SP^N^)  by  the  permutations  of 

seems  to  be  complicated  for  general  k  ;  the  difficulties  are  illustrated 

below  in  the  case  k  =  3  0 

Suppose  a  permutation  jr  of  N  has  j„  cycles  of  length 

m  l 

i  and  let  JL  denote  the  set  containing  the  ij^  nodes  involved  in 
these  cycles  for  i  =  l,2,,.0,m  „  We  begin  by  determining  the  edge 
cycles  induced  in  SP^(j_.  )  by  it  . 

Case  Is  If  i  =  0  (mod  3)  ♦  the  edge  cycles  of  SP„(J\)  are  of 

' — — ' —  j  * 

two  kinds;  first,  there  are  edge  cycles  which  involve  only  edges  having 

i/3 *  2i/3 , 

the  form  {a,  jt  n  '  ^(a)}  ,  a  c  J,  ,  and  second,  edge  cycles 

which  involve  only  edges  having  the  form  (a,b,c)  ,  a,b,c  €  ,  not 

of  the  first  kind,  There  are  ij  /3  edges  of  the  first  kind  and  these 
edges  split  into  cycles  of  length  i/3  ,  so  there  are  j.  edge  cycles 
of  length  i/3  induced  by  st  among  these  edges.  The  remaining  edges 

es  edge  cycles  of  length  i 
-  ij^/3j"  edge  cycles  of 

length  i  . 


are 


indui 


ij./3  in  number  and  jt  i 


among  these  edges,  so  there  are 


Case  2:  If  i  =  1,2  (mod  3)  the 

into  cycles  of  length  i  ;  hence, 
length  i  induced  in  SP^(J^)  by 


edges  in  SP^(J_^)  split 


there  are 


i/lJi 


edge  cycles  of 


jt 


'  • 

13.  V^vla  2S,' 
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three  cases  depending  on  the  parity  of  i^  and 


Case  3»  If  i-^  and  i  are  both  odd,  then  the  edges  having  the  form 


{a,b,c}  with  a  e  J.  ,  b,c  e  J.  ,  are  i.j. 

Li  x2  1 


Vi. 


in 


number,  and  they  split  into  cycles  of  length  ^east  coinmon 

multiple  of  i^  and  ±2  .  Thus,  there  are  $±  j£  (i1,i2)(i2ji  ~1)/2 

edge  cycles  of  this  type$  interchanging  the  subscripts  1  and  2  gives 

the  number  of  remaining  edge  cycles  also  of  length  [i  ,1^]  .  Altogether 

there  are  j  j  (i.,i  )(i  j  +  i  j  -2)/2  edge  cycles  of  length 
12  1  d  2 

[i.,i  ]  induced  by  jt  among  the  edges  of  SP  (j.  U  J  )  -  SP  (J.  )  - 

12  j  lg  j  1^ 

SP  (j.  )  when  i1  and  ip  are  both  odd, 

j  1  c- 


Case  h:  If  i,  and  i^  are  of  opposite  parity,  say  i0  is  even,  the 

—————  12  2 

edge  cycles  induced  by  jt  are  of  three  kinds£  first,  those  involving 

V2 

edges  having  the  form  {a,b,  jc  (b ) }  with  a  e  J.  ,  b  e  J.  ’ 

12 

second,  those  involving  the  edges  having  the  form  {a,b,c}  with  a  €  J.  , 

V2  1 

b,c,  e  J.  ,  c  ^  «  (b)  ;  third,  those  involving  edges  having  the  form 

L2 


[a,b,c]  with  a,b  e  J,  ,  c  e 


The 


re  are  i-j.  °  ip  j  .  /2  edges 

1  ‘2  1  2 

of  the  first  form  and  these  split  into  cycles  of  length  [i^,ig/2] 

hence,  there  are  j.  j.  .  (i  ,i  /2)  edge  cycles  of  length  [i  ,i  /2] 

11  L2  1  2 

induced  in  this  set.  The  edges  of  the  second  form  are  i^j^  fgj^  ^2^i 
in  number  and  they  form  cycles  of  length  [i^,ig]  hence,  there  are 


2  )/2 


' 


. 


it  >  o  , 


■ 
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j.  j.  (i-.  >  ip)  (ipj .  -2)/2  cycles  of  length  [1  -  ,  i  ]  induced  in  this 
12  2  1  2 


set.  Finally,  the  ip j . 

12 


i,  J  • 
lJi. 


edges  having  the  third  form  split 


into  cycles  of  length  [  i  i  ]  as  well,  so  there  are  j.  j.  (i  ,i  )(i  j.  -2)/2 

■*-  ig  l  £_  1  i^  ' 

cycles  of  this  length  induced  among  the  edges  in  the  last  set. 


Summing  up  these  results,  we  have  that  if  i^  is  odd  and 

ip  is  even,  then  there  are  j  j  (i  ,i  )/2  edge  cycles  of  length 

1  2'  7 

[il^^/2]  and  (il,i2)(il^i  +  ^2^i  “^V2  cycles  of  length 

[i  ,i  ]  induced  in  the  set  SP  (j.  U  J.  )  -  SP  (j.  )  -  SP  (j.  ) 

^  L1  X2  ^  11  ^  12 

by  it  . 


Case 


Jj;  If  i^  and  ip  are  both  even,  we  can  ignore  the  parity  of 


i^  for  the  moment  and  determine  the  edge  cycles  among  the  edges  having 

ip/2 

the  forms  {a,b,jt  (b)j,  a  e  J.  ,  b  e  J.  ,  and  {a,b,cj,  a  e  J.  ,  b,c,  e  J.  , 

11  12  11  12 

ip/2  ^ 

c  /  jt  (b),  just  as  in  Case  4  above.  This  gives  2  j  „  j.  (i..,ip/2)  edge 

1112 

cycles  of  length  [i^i^/2]  and  ji  ji  (i^  ig)  (i2,  ^  -2)/2  edge  cycles 

of  length  [i^i^].  By  interchanging  the  subscripts  1  and  2  we 
obtain  the  number  of  remaining  edge  cycles.  In  conclusion,  if  i  and 
i  are  both  even,  then  there  are  j.  j„  (i1,ip/2)  cycles  of  length 

C.  1  Ip  t— 

[ ip  ip/2]  >  JjL  (ij_/2,  i2 )  cycles  of  length  [i^/2,  ig]  ,  and 

ji  jt  (i1,ip)(i1ji  +  ip j i  -4)/2  cycles  of  length  [i^ip]  induced  by 

tc  among  the  edges  of  SP  (J.  U  J.  )  -  SP  (j.  )  -  SP  (j,  )  , 

3  11  l2  ^  11  ^  l2 
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Finally,  there  are  edge  cycles  induced  by  n  involving  nodes 
chosen  from  three  different  sets  J\,  . 

Case  6:  There  are  i  j  i  j.  i  j,  edges  having  the  form  {a,b,c} 
with  a  e  J.  ,  beJ.  ,  ceJ.  and  jt  induces  i. j .  ip j .  i,j. 

il  i2  11 1  ^  L2  ^  3 

/[il,i2,i^]  cycles  of  length  [i^,ig,i  ]  among  these  edges. 

^1  ^2  ^1  ^2 
If  we  write  c(jt)  =  f^  £*  ...  and  c(jt*)  =  g1  g2  , 

where  jt  is  a  permutation  of  and  jt*  is  the  permutation  induced 

by  jt  in  SP  (Nm)  ,  then  we  have  the  following  correspondence  between 

the  factors  of  c(jt)  and  c(jt*)  2 
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When  applying  these  formulas,  the  pair  or  triple  of  sub¬ 
scripts  i.,i0  or  i  ,i  , i  is  considered  exactly  once.  Using 
1  d  i  d  3 

Riordan's  [19]  list  of  cycle  indices  for  the  symmetric  groups  on 
Nm  and  the  correspondences  given  in  (ll),  (12),  and  ( 13) ,  we  get 
the  following  cycle  indices  for  the  groups  induced  in  SP^(Nm)  for 
m  =  6,7,8,  and  9  >  respectively. 


1  r  20  1C  8  6  on  2  6  2  2  2  ic  10 

(14)  {§!  +  158^2  +  ^5glg2  +  80glg3  +  120^^S  +  15g2 

+  lSOgggJ  +  120S2g6  +  I44g^  } 

(15)  JT  {g^5  +  21g^5g20  +  105g^gg^  +  70g^g^°  +  420g^g2g^ 

+  105g^26  +  280g^g^  +  840g1g^gJ  +  210g1g2g^ 
+  420g1g1^g6g^2  +  840g2g^  +  504g^  +  504g^g2Q 
+  720g^  ) 


+  28g^6g25  +  210g^2g22  +  l^g^gj5  +  420g^g25 

s, 

+  1120g^g|gjgg  +  420g^gggJ1  +  l680g^gggjgg 

+  1120g^g®g|  +  1120g^gj8  +  2520g^g|g“ 
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+  +  2688g1g^g^5  +  i*032g1g^g^0 

+  I344g1gj1  +  105g28  +  ^Og^gJ1  +  3360g2g9 

+  3360g2g2g8  +  126ogJ4  +  5?6og®  +  501+Ogg  ) 


(17) 


(sf  +  SSgfgl1  +  378gfgf  ♦  I68g21  sf  +  I260g;°gf 
+  2520g9g^gl1g5  +  756g}°gJ6g|  + 

+  7560g^g2g^g®  +  3360g^g|7  +  7560gJgggJ6 
+  3024gig^6  +  2520g]g2°g5g-?  +  10080g^gj5g^ 

+  11340g2g9gJ6  +  15120g1ggg^gj|g6gi2  +  l8l^ltg^g2g®g^0 

+  lOOSOgjSggjgg2  +  2280g7g27  +  15120glg2gjg|jgggi2 

p  k  f,  p  pn  k  k 

+  9072g2g5g10  +  11340g2g4  +  24l92g1g5g^gJ'5 

+  3021+Og^ggjgg5  +  25920g^2  +  lSlltltg^g^g^g^ 

+  25920g°8^  +  45360g^gg°  +  40320g^2  ) 


The  number  of  non- isomorphic  (3,S)-graphs  with  n  edges  and 

m  nodes,  for  m  =  6, 7, 8,9*  and  n  =  0,1, ...,m  ,  is  the  coefficient  of 

xn  in  (14),  (15),  (l6),  and  (17)  respectively,  when  1  +  x"*'  is 

substituted  for  g.  in  each  of  these  relations. 
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